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1 INTRODUCTION 3 

1 Introduction 

One of the outstanding problems in particle physics is the determination of the 
structure of hadrons such as the proton and neutron in terms of their fundamental 
quark and gluon degrees of freedom. Over the past twenty years two fundamentally 
different pictures have developed. One, the constituent quark model is closely 
related to experimental observation. The other, quantum chromodynamics is based 
on a covariant non-abelian quantum field theory. The front form (also known as 
light-cone quantization) appears to be the only hope of reconciling these two. This 
elegant approach to quantum field theory is a Hamiltonian gauge-fixed formulation 
that avoids many of the most difficult problems in the conventional equal-time 
formulation of the theory. 

The natural gauge for light-cone Hamiltonian theories is the light-cone gauge 
A"*" = 0. In this physical gauge the gluons have only the two physical transverse 
degrees of freedom. One imagines that there is an expansion in multi-particle oc- 
cupation number Fock states. But even in the case of the simpler abelian quantum 
theory of electrodynamics very little is known about the nature of the bound state 
solutions in the strong-coupling domain. In the non-abelian quantum theory of 
chromodynamics a calculation of bound-state structure has to deal with many dif- 
ficult aspects simultaneously. Confinement, vacuum structure and chiral symmetry 
inter-twine with the difficulties of describing a (relativistic) many-body system and 
the non-perturbative renormalization of a Hamiltonian. 

In the conventional approach based on equal-time quantization the Fock state 
expansion becomes quickly intractable because of the complexity of the vacuum. 
Furthermore, boosting such a wavefunction from the hadron's rest frame to a 
moving frame is as complex a problem as solving the bound state problem itself. 
The presence of the square root operator in the equal-time Hamiltonian approach 
presents severe mathematical difficulties. 

Fortunately 'light-cone quantization' offers an elegant avenue of escape. It can 
be formulated independent of the Lorentz frame. The square root operator does 
not appear, and the vacuum structure is relatively simple. There is no spontaneous 
creation of massive fermions in the light-cone quantized vacuum. 

In fact, there are many reasons to quantize relativistic field theories at fixed 
light-cone time. Dirac |]I| showed, in 1949, that in this so called 'front form' of 
Hamiltonian dynamics a maximum number of Poincare generators become inde- 
pendent of the interaction, including certain Lorentz boosts. In fact, unlike the 
traditional equal-time Hamiltonian formalism, quantization on a plane tangential 
to the light-cone, on the 'null plane', can be formulated without reference to a 
specific frame. One can construct an operator whose eigenvalues are the invari- 
ant masses of the composite physical particles. The eigenvectors describe bound 
states of arbitrary four-momentum and invariant mass, and allow the computation 
of scattering amplitudes and other dynamical quantities. In many field theories 
the vacuum state of the free Hamiltonian is also an eigenstate of the light-cone 
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Hamiltonian. The Fock expansion built on this vacuum state provides a complete 
relativistic many-particle basis for diagonalizing the full theory. 

The main thrust of these lectures will be to discuss the complexities that are 
unique to this formulation of QCD, in varying degrees of detail. The goal is to 
present a self-consistent framework rather than trying to cover the subject ex- 
haustively. A review all of the successes or applications will, however, not be 
undertaken. 

One of the reasons is, that the subject was reviewed recently 0. Another is 
that other lecturers in this school emphasize complementary aspects. Stan Brod- 
sky shows how the knowledge of the the light-cone wavefunction has impact on 
hadronic physics and exclusive processes. Steve Pinsky demonstrates how the 
method of discretized light-cone quantization is constructive for analyzing super- 
symmetric string and M(atrix) theories. Last not least, Simon Dalley expands on 
the transverse-lattice calculations within the light-cone approach. 

Comparatively little space will be devoted to canonical field theory, just so much 
as to plausibilize that a light-cone Hamiltonian exists. This so called 'naive Ha- 
miltonian' will be derived and written down explicitly as a Fock-space operator in 
the light-cone gauge and disregarding zero modes. For historical and paedagogical 
reasons these notes will be rather outspoken for one space and one time dimension, 
mostly to show that periodic boundary conditions (discretized light-cone quanti- 
zation) are helpful, indeed, for solving the bound-state problem in a relativistic 
theory. The attempt is made to be complementary to the review and some new 
material on the Schwinger model is included. 

The bulk of these notes deals with the many-body aspects of a gauge field theory 
in the real world of 3+1 dimensions. A hadron not only contains the valence quarks 
but also an infinite amount of gluons and sea-quarks. Progress often comes with 
new technologies. The presentation of the method of iterated resolvents therefore 
takes broad room. It allows to derive a well-defined effective interaction. Some thus 
far unpublished work is included, in particular more instructive examples. The 
presentation is separated into two parts. In the first considerations are essentially 
exact. In the second some approximations and simplifications are admitted and well 
marked in the text. One arrives at the effective interaction in the form of a tractable 
and solvable integral equation. Its solutions allow to construct explicitly the many- 
body amplitudes corresponding to see-quarks and gluons by comparatively simple 
quadratures, without the need of solving another bound-state problem. Explicit 
equations for that are given. Moreover, some new research work will be included 
in the section on renormalization. 

What is not included in these notes, however, is a complete survey of the lit- 
erature. Mostly for the reasons of space, it is refered to the some 469 items of 
0. I apologize with my colleagues whose work is not mentioned. But I will be 
careful to cite the work which I need for the present discussion and presentation. 
Some selected monographies which I found useful to consult are [^, ^, ^ ^, ^ , some 
selected review articles or conference proceedings might be H, 0, [TO, |ll[ . 
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2 Forms of Hamiltonian Dynamics 

Dirac defines the Hamiltonian if of a closed system as that operator whose action 
on the state vector | t ) has the same effect as taking the partial derivative with 
respect to time t, i.e. 

H\t)=z^^\t). (1) 

The concept of an Hamiltonian is applicable irrespective of whether one deals 
with the motion of a non-relativistic particle in classical mechanics or with a non- 
relativistic wave function in the Schrodinger equation, and it generalizes almost 
unchanged to a relativistic and covariant field theory. 

In a covariant theory, the very notion of 'time' becomes, however, questionable 
since the time is only one component of four- dimensional space-time. But the 
concept of usual space and of usual time can be generalized. One can define 'space' 
as that hypersphere in four-space on which one chooses the initial conditions. The 
remaining fourth coordinate can be understood as 'time'. 

More formally, one conveniently introduces generalized coordinates x'^. Starting 
from a baseline parametrization of space-time like the instant form in Figure |l], one 
parametrizes space-time by a coordinate transformation x" = x'^{x'^). The metric 
tensors for the two parametrizations are then related by 

The physical content of the theory can not depend on such a re-parametrization. 

But the raising and the lowering of Lorentz indices are then non-trivial opera- 
tions. As an example consider the front form parametrization in Figure |l|. Interpret 
x^ = x^ — x^ as the new space coordinate and denote it by x'^ = ct — z. Then 
go _ 2;0_|_a;3 ]2iust be interpreted as the new time coordinate denoted by x~^ = ct+z, 
or by the 'light-cone time' t = t + z/c. Of course, one also could exchange the 
two. Since the lowering operation is x^ = g^uX" , both x+ = |x~ are space, and 
x_ = ^x~^ are time coordinates. The new space derivative is therefore d- = ^d~^, 
while (9+ = |9~ is a time-derivative. The Lorentz indices '+' and '-' have a different 
physical meaning, depending on whether they occur up- or down-stairs. Co- variant 
and contra- variant vectors are different objects. The Hamiltonian is only one com- 
ponent of a four- vector P^, particularly its time-like component. Taking the partial 
time derivative of the state vector 

defines then P+ = |P^ as the Hamiltonian in the transformed coordinates, in line 
with Eq.(^, and P = ^P^ as the longitudinal momentum. 

Following Dirac |]T| there are no more than three different parametrisations of 
space-time. They are illustrated in Figure |l], and cannot be mapped onto each 
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Figure 1: Dirac's three forms of Hamiltonian dynamics. 



other by a Lorentz transform. They differ by the hypersphere on which the fields 
are initiahzed. They have thus different 'times' and different 'Hamiltonians'. Dirac 
|1[] speaks of the three forms of Hamiltonian dynamics: The instant form is the 
famihar one, with its hypersphere given by t = 0. In the front form the hypersphere 
is a tangent plane to the light cone. In the point form the time-like coordinate is 
identified with the eigentime of a physical system and the hypersphere has a shape 
of a hyperboloid. 

Which of the three forms should be prefered, is an ill-posed question, since all 
three forms must yield the same physical results. Comparatively little work has 
been done in the point form. The bulk of research on field theory implicitly uses 
the instant form. Although it is the conventional choice for quantizing field theory, 
it has many practical disadvantages. For example, given the wavefunctions of an 
n-electron atom at an initial time t = 0, ipni^iii = 0), one can use the Hamiltonian 
H to evolve ipnixijt) to later times t. However, an experiment which specifies the 
initial wave function would require the simultaneous measurement of the positions 
of all of the bounded electrons. In contrast, determining the initial wave function 
at fixed light-cone time r = only requires an experiment which scatters one plane- 
wave laser beam, since the signal reaching each of the n electrons, along the light 
front, at the same light-cone time r = tj + Zi/c. 

Dirac's legacy had been forgotten and was re-invented several times. The front 
form approach carries therefore names as different as Infinite- Momentum Frame, 
Null-Plane Quantization, Light- Cone Quantization, or most recently Light-Front 
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Quantization. In the essence they are all the same. The infinite-momentum frame 
is a misnomer since the total momentum is finite and since the front form is frame- 
independent and covariant. Light-cone quantization is also unfortunate, since the 
inital data are set one a plane tangential to but not on the light cone, and since 
both equal-time or equal light-cone-time quantization stand both for the same 
quantization, for a quantum as opposed to a classical theory. We propose to stay 
with Dirac's different 'forms of Hamiltonians'. 

2.1 The canonical Hamiltonian for gauge theory 

The prototype of a field theory is Faraday's and Maxwell's electrodynamics. Every 
field theory has its own canonical Hamiltonian and is governed by the action density. 
The Lagrangian is the subject of the canonical calculus of variation, given in many 
text books 0. Its essentials shall be recalled briefiy. 

The Lagrangian, in general, is a function of a finite number of fields (pri^) 
and their first space-time derivatives dfj,(f)r{x), thus C = C [(f)r,dfj,(f)r]- Independent 
variation of C with respect to (pr and (9^0,- results in the equations of motion, 

A/* 

d^TT^ - 6C/6(l)r = , with <[0] 



6 {d^(t)r) ■ 



canonically refered to as the Euler equations. The canonical formalism is partic- 
ularly suited for discussing the symmetries of a field theory. Every continuous 
symmetry of the Lagrangian is associated with a vanishing four-divergence of a 
current and a conserved charge. Since C does not explicitely depend on the coor- 
dinates, every field theory in 3-1-1 dimensions has ten conserved four-currents. The 
four-divergences of the energy-momentum tensor T^^ and of the boost-angular— 
momentum stress tensor J-^'^*^ vanish, 

dxT^" = 0, dxJ^'^" = 0. 

As a consequence the Lorentz group has ten 'conserved charges', the 4 components 
total momentum P'^ and the 6 components of boost-angular momentum M^^", 

P" = / dtooT'"' , and M"" = f dcuoJ^'"". (2) 

Jn Jn 

With the totally antisymmetric tensor in four dimensions ex^up (^0123 = 1), the 
three-dimensional surface elements of a hypersphere are dux = ex^.updx^dx'^dx^/Sl. 
A finite volume is thus fl = J duj^ = J dx^dx'^dx^, and correspondingly in the front 
form du^ = J dx+d'^x±. The time-like components of P^ is the Hamiltonian, i.e. 
Pq for the instant form and P+ for the front form. The transition from instant to 
front form is thus simple: substitute '0' by '+'. 

Working out the canonical procedure for QED with its Lagrangian, 

C = -^F^^F^, + ^ [^ {lYD^ -m)^+ h.c. ] , (3) 
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where F^^" is the electro-magnetic field tensor and D^ the covariant derivative, 



F"^ = d^A"" - d^A", 



D, 



da - igA 



^l1 



one ends up with a manifestly gauge-invariant total momentum, 

Jiu^ [f'^F^, + ^g'.F^'F,^ + 1 [2^7°/},^ + h.c. 



dujA 



F+^F^, + i^+F'^^F^A + ^[2^7+^.^ + h-c. 



in both the instant and the front form, respectively. The boost angular momenta 
will not be used explicitly. 

The gauge invariant Lagrangian density for QCD is 

C = -iTr(F^'^F^,) + ^[^(z7'^D^-m)^+ h.c], 

= -^FrF;, + i[^(z7^D,-m)vl/+ h.c]. 

The color-electro-magnetic fields and the covariant derivative are now 

F'^^ = 9^A'^-9''A^ + i^[A^,A''] = T" (9M^ - 9^A^ - ^r"M(f A^) , 



D^,, = 5,,d^ + igA''^, 



<5ec'9^ + igT:,A^^. 



As compared to QED, each local gauge field A^^{x) is replaced by the 3x3 matrix 
A^(a;). All such matrices can be parametrized A^' = T^^iA'^. More generally for 
SU(N), the vector potentials A^ are hermitian and traceless N x N matrices. The 
color index c (or <^) runs now from 1 to tt-c, and correspondingly the gluon index 
a (or r,s,t) from 1 to ra^ — 1. No distinction will be made between raising or 
lowering them. In order to make sense of expressions like \1'A'^\1' the quark fields 
"^caix) must carry a color index c. They, as well as the Dirac indices, are usually 
suppressed. The color matrices T"^/ obey 



yr y. 



i/"'Xc' and Tr {T^'T') = -6'^. 



(4) 



For SU(2) the color matrices are T"" = \(J°'-, with a"' being the Pauli matrices. 
The structure constants /'''^* are therefore the totally antisymmetric tensor erst- 
For SU(3), T" = ^A** with the Gell-Mann matrices A*^, with the corresponding 
structure constants tabulated in the literature P]. Everything proceeds in analogy 
with QED. The energy-momentum vector, 

^du, {F'^^F:, + \glF:^F:^ + \\i^^'T'^Dl^ + h.c 



Pu 



dujA 



K^Fi + \gtF:^Kx + \ [i^i'-T'-Dy + h.c. 



(5) 



is manifestly gauge-invariant both in the instant and the front form. 
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2.2 The Poincare symmetries in the front form 

The ten constants of motion P^ and iVf"^ are observables with real eigenvalues. It 
is advantageous to construct representations in which the constants of motion are 
diagonal. But one cannot diagonalize all ten constants of motion simultaneously 
because they do not commute. The algebra of the four-energy-momentum P'^ = p^ 
and four-angular-momentum M^'^ = x^p'^ — x'^p^ for free particles with the basic 



commutator j^[x'^,pi,] = 6^ is 




ih 


= gf^P" - g'^'P" , 4 \P^, P^ 

in 


and ^\MP%Mn = 
in 


= gf'M'''' + g^^Mf" - gPf'M'"' 



It is postulated that the generalized momentum operators satisfy the same com- 
mutator relations as a single particle. They form thus a group, the Poincare group. 

It is convenient to discuss the structure of the Poincare group in terms of the 
Pauli-Lubansky vector V^ = e'^^^'^P\M^y. V is orthogonal to the generalized mo- 
menta, P^V^ = 0. The two group invariants are the operator for the invariant 
mass-squared M^ = P^P^ and the operator for intrinsic spin-squared V"^ = V^V^. 
They are Lorentz scalars and commute with all generators P^ and M^^ , as well 
as with all V^. A convenient choice of six mutually commuting operators is: The 
invariant mass squared M^ = P'^P^, the three space-like momenta P~^ and P±, the 
total spin squared S"^ = V^V^, and one component of V , say V^ = S^- 

Inspecting the definition of boost-angular momentum M^^, in Eq.(H) one identi- 
fies which components are dependent on the interaction and which are not. Dirac 
|T|] calls them complicated and simple, or dynamic and kinematic, or Hamiltonians 
and Momenta, respectively. In the instant form, the three components of the boost 
vector Ki = Mjo are dynamic, and the three components of angular momentum 
Ji = ^ijkMjk are kinematic. As noted already by Dirac, the front form is special in 
having four kinematic components of M^y {M^_, M12, Mi_, M2-) and only two dy- 
namic ones (M+i and M+2)- In the front form one deals thus with seven mutually 
commuting operators 

(M+_, M12, Mi_, M2-), and all P^, 

instead of the six in the instant form. These symmetries imply the very important 
aspect of the front form that both the Hamiltonian and all amplitudes obtained 
in light-cone perturbation theory are manifestly invariant under a large class of 
Lorentz transformations: 



p+ ^ C|| p+ , p_L ^ R , P ^ C*!! V 



■y2 



p+ — > p+ , p± -^ Pl+ p'^C± , p -^ p +2pi ■ C±+ p'^Cl , 

-I- + ->2 -^2 

p+ ^p+ ^ a: -^ p± , 

i.e. for parallel boosts, transverse boosts, and rotations, respectively. All of these 
hold for every single particle momentum p^, and for any set of dimensionless c- 
numbers Cm and Ci . 
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2.3 The light-cone Hamiltonian operator 

The four- vector of energy- momentuni for gauge theory in Eq.(^) contains time- 
derivatives and other constraint field components. They will be eliminated in this 
section using the equations of motion with the goal to express P" in terms of the 
free fields and to isolate the dependence on the coupling constant. 

The color-Maxwell equations are four equations for determining the four func- 
tions y4^. One of the equations of motion is identically fuUfilled by choosing the 
light-cone gauge A'^ = 0. Two of the equations give expressions the time deriva- 
tives of the two transversal components A1. The fourth is the Gauss' law in the 
front form, df^F^~^ = gJ^, or explicitly 

- d^d^A- - 9+9,4, = gJ^- (6) 

It contains only space-derivatives and is a constrained equation for the components 
of A^. For the free case {g = 0), A^ reduces to A^ and therefore to 



^a — y^a 5 ^-La5 -^a) — I 0' ^^a.-, ~ a~ 



diA\ 



As a consequence, A^ is purely transverse. The formal inversion of Eq. (|^) is there- 
fore 

which must be substituted everywhere. The inverse space derivatives (z9"'")~ and 
[id^Y , used here and below, are actually Green's functions. Since they depend 
only on a;~, they are comparatively simple. 

The color-Dirac equations (z7^D^ — m)^ = can be used to express the time 
derivatives 9+\E' as function of the other fields. After multiplication with 7° one 
gets first 

(z7Vr"^+ + i7%~r"D^ + ic^J'^Dl^^ = rnp^H, 

with the usual Dirac matrices /3 = 7° and a'' = 7*^7*''. In order to isolate the time 
derivative one introduces the projectors A-t and projected spinors \I'± by 

A± = -(l±«^) and ^± = A±^ . 
Multiplying the color-Dirac equation once with A+ and once with A_ gives 

and 2i9_^_ = (m/?-iaiT"D^i)^+ + 2^A^T"^_, (8) 

a coupled set of spinor equations. The first of them contains a time derivative. The 
second contains a space derivative and is a constraint equation. The component 

^- = ^ {m(3 - lalT^Dl) ^^ (9) 
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must therefore be substituted everywhere. The time derivative becomes then 
2z9+^+ = 2gAlT^^+ + (m(3 - laiT^Dl^ -^ (m/3 - ialT^Dl^ ^+ . 
In analogy to the color-Maxwell case one defines free spinors by 

^ ^ 2Z(7_ 

Contrary to the full spinor, ^ is independent of the interaction. 

Inserting the above expressions into Eq.(^, the space-like components of P^ 
become 

Pk = jdx+(fx^ (f 7""^^^^ + K -9+9^1^) , for fc = 1, 2, -. (10) 

Inserting them into P+ gives rather lengthy expressions, which are conveniently 
written as a sum of five terms 

P+ = T + V + Wi + W2 + W^. (11) 

Only the first term survives the limit (7 ^ 0, and therefore is called the free part 
of the Hamiltonian, or its 'kinetic energy' 



T = \ \dx^dW (^7-^ ""' ^ir^^' ^ + ~K^^^fK) ■ 



The vertex interaction 



V = g Jdx^d'x^ r,Al, with j:{x) = ^^''T'^^ + r'^d^AlAl, (12) 

is linear in the coupling constant and is the light-cone analogue of the conventional 
J^A'^-structures in the instant form. Note that the current J^ has contributions 
from both quarks and gluons. The four-point gluon interaction 

W, = ^^ jdx^dW B^,''B%, with i?r = /"^^Ir^r, 



describes the four-point gluon- vertices which is quadratic in g. The remainders are 
the 'instantaneous interactions'. They are characterized by the inverse derivatives. 
The instantaneous gluon interaction arises from the Gauss equation, 

W2 = — fdx+d^x^ J+—^J^, 

and is the light-cone analogue of the Coulomb energy. The instantaneous fermion 
interaction originates from the light-cone specific decomposition of Dirac's equation 

W, = ^Jdx+d'x^ ^YT^A; ^ (YT'Al^) . 
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It has no analogue in the instant form. 

Most remarkable, however, is that the relativistic Hamiltonian is additive in 
the 'kinetic' and the 'potential' energy, very much like a non-relativistic Hamilto- 
nian H = T + U. In this respect the front form is distinctly different from the 
conventional instant form. 

2.4 The free field solutions 

The free solutions of the Dirac and the Maxwell equations are in the front form 






dp'^cPp± 
A^ ^2p+{27r) 



The properties of the Dirac spinors Ua and Va, and of the polarization vectors e^, 
are given for example in . The single particle state are specified by the quantum 
numbers q = {p'^, p±x, P±y, A, c, /). Their creation and destruction operators are 
subject to the usual relations 

a(g),at(g')] = {biq),b\q')} = {d{q), d\q')} = 5{p^ - p+ ')5^'\p^ - p[_)5^; St 5^; , 

which carry the operator structure of the theory. When inserting the free fields into 
P/i one can integrate out the dependence on x^, producing essentially Dirac delta- 
functions in the single particle momenta, which reflect momentum conservation: 







In detail this can be quite laborious, as shown by the example with the fermionic 
contribution to the vertex interaction 



Vf = g jdx+d^x^ m{x)-i^T''^!{x)Al{x 



a;+=0 



g fdpjdy±i fdp2d^p±2 fdpjd^p±3 



X 



y(27r)3 Ai,A2,A3Ci,C2,a3-^ ^2pt "^ ^J2p^ "^ ^J2p^ 

dx^dW r^^t(^^)^^(p^^ Ai)e+^^^^' + d{q,)v^{p,, Ai)e-*P^^') T^^ 



(27r)3 

X 7^;3 {d\q,)vp{p2, A2)e+^^^^ + b{q2)up{p,, A2)e-^P^^) 

X (a^(g3)e;(P3, A3)e+*^«-' + a(g3)e^(P3, A3)e 



Cl,C2 



-vp-iX 



Note that the sum of these single particle momenta is essentially the sum of the 
particle momenta minus the sum of the hole momenta. Consequently, if a particular 
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V4 



+ 



+ 



A 



V 



"-1 "-2 "-3 

Ay 

i.+ i.+ i.+ 
"-1 "-2 "-3 



0203 



A 



+ 



"'I "'2 "-3 

i,+ i.+ i.+ 

Aj-| AjO '»jQ 

i.+ i.+ i.+ 

At]^ "-2 "^3 






Table 1: The vertex interaction in terms of Dirac spinors. The matrix elements 
Vn are displayed on the right, the corresponding (energy) graphs on the left. All 
matrix elements are proportional to Ay = gS{ki\k2+3)S^'^\k±^i\k±^2 + ^±,3); with 
g = g/J2{2Tr)^. For the periodic boundary conditions one uses g = g/^/Ti. 



term has only creation or only destruction operators as in 



its contribution vanishes since the light-cone longitudinal momenta p^ are all posi- 
tive and can not add to zero. As a consequence, all energy diagrams which generate 
the vacuum fluctuations in the usual formulation of quantum field theory are absent 
in the front form. 



2.5 The Hamiltonian as a Fock-space operator 

The kinetic energy T becomes a sum of 3 diagonal operators 



T = fdk.d% Y. ( 

A,c,/ 

-m"^ + kf 



E 



A;_ 



m^ + kl 
k_ ^ g 



h\hq + d\dq + a\aq 



Here and in the sequel it is convenient to abbreviate the integration and summation 
over the single particle coordinates by the symbol Yl, to replace for instance h{q) 



with hq. 



The vertex interaction V becomes a sum of 4 operators 

V = V1 + V2 + VS + V4 
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^ 



^ 



i-j>y-^ji-j-^ 



iK/CC'v' 2 






— <J 



-^3,1 
-^3,2 

-^5,2 
-^5,3 
-^6,1 
-^6,2 



2A (miT'^7+M2) (t^37+T"M4) 



"'1 "-2 "-3 

A 



Ra 



(MiT'^V47+/3T"'M2) 



,+A;.+fctfc+ (^f - kt) 



2ktA 



3 '^4 



(iiir'^7+M2) (estC^i) 



"'1 "-2 "-3 



A 



(kt - ktf 

(t;3T"Vt7+/2T'^^M4) 



+ i,+h+h+ (K - kt) 



"-1 "-2 "-3 

A 



(vsT'^^hl^fiT'^'Ui) 



tktktkt (^1 - ^4 



At| At2 At3 
^1 /Ci ~r 1^2 



)A(i;3T"7+M4) (el^C^es) 



ktk+k+kt (^1 - ^2+)' 



Hi ^2 /t3 

2A 



(A;+ - fc2+)2 



i.+ J,+ J,+ ^H 
"1 "-2 "-3 "-4 



(ete3) (e2e4) C:^a,C:, 



a-ja^ 



Table 2: The fork interaction in terms of Dirac spinors. The matrix elements Fnj 
are displayed on the right, the corresponding (energy) graphs on the left. All matrix 
elements are proportional to A = g'^5{kf\k2 + k^ + k^)5^'^\k±^i\k±^2 + ^±,3 + ^±,4)7 
with g = g/ J2{2tx)^. For the periodic boundary conditions one uses g = g/yVt. 



= E [^i^sas Vi{l- 2, 3) + h.c] + E [4^203 1^2(1; 2, 3) + h.c." 
1,2,3 1,2,3 

+ E W2h V^3(l; 2, 3) + h.c] + E H«2a3 ^4(1; 2, 3) + h.c. 



+ 



1,2,3 



1,2,3 



It connects Fock states whose particle number differs by 1. The matrix elements 
V^(l; 2, 3) = Vn{qi; q2, (I3) are simple functions of the three single-particle states g^, 
which are given in Table ^ 

The four-point interactions are broken up conveniently into fork and seagull 
interactions, F and S, depending on whether they have an odd or an even number 
of creation operators, thus 

P^ = T + V + F + S. 

The fork interaction F becomes then a sum of 6 operators, 

F = F1 + F2 + F3 + F4 + F5 + F6 
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= J2 [&ife24fe4i^i(l;2,3,4) + h.c.] + [4c?2&34i^2(l;2,3,4) + h.c." 

1,2,3,4 

+ Y. [^1^20304 Fsil; 2, 3, 4) + h.c] + [40^20304 ^4(1; 2, 3, 4) + h.c. 
1,2,3,4 

+ Y. [aWdsh F^il; 2, 3, 4) + h.c] + [4030304 i^6(l; 2, 3, 4) + h.c] .(13) 
1,2,3,4 

It changes the particle number by 2. The matrix elements are given in Table Q. 
The seagull interaction S, finally, becomes a sum of 7 operators 

= E b\blhb,S,{l,2;3,4)+ E 44^^3^^4 52(1, 2; 3, 4) 
1,2,3,4 1,2,3,4 

+ E &14M4 53(1,2;3,4)+ E &I4M4 54(1,2;3,4) 
1,2,3,4 1,2,3,4 

+ E 444^4 55(1, 2; 3, 4)+ E (^i4a3a4 56(l,2;3,4) + h.c) 
1,2,3,4 1,2,3,4 

+ E 44«3a4 5'7(l,2;3,4). 
1,2,3,4 

Its matrix elements are given in Table |^. It can act only between Fock states with 
the same particle number. 

The above results are are quite generally applicable: They hold for arbitrary 
non-abelian gauge theory SU{N). They hold for abelian gauge theory (QED), 
formally by replacing the color-matrices T^^, with the unit matrix and by setting 
to zero the structure constants f"'^'^, thus B'^" = and x^ = 0. They hold for 1 
time dimension and arbitrary d + 1 space dimensions, with i = 1, . . . ,d. All what 
has to be adjusted is the volume integral Jdx+dx±^i . . . dx±^d- 

3 The hadronic bound-state problem 

One has to find a language in which one can represent hadrons in terms of relativis- 
tic confined quarks and gluons. As reviewed in |T^, the Bethe-Salpeter formalism 



has been the central method for analyzing hydrogenic atoms in QED and provides 
a completely covariant procedure for obtaining bound state solutions. However, 
calculations using this method are extremely complex and appear to be intractable 
much beyond the ladder approximation. It also appears impractical to extend this 
method to systems with more than a few constituent particles. 

An intuitive approach for solving relativistic bound-state problems would be to 
solve the instant form Hamiltonian eigenvalue problem 



for the hadron's mass and wave function. Here, one imagines that |\E') is an expan- 
sion in multi-particle occupation number Fock states, and that the operators H 
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1 -• 5 —3 



A (miT"7+M3) (M27+y"^4) 

hj^ rv2 "-'3 'i-'4 ^ ^ ^ ^ 

2A (miT'*7+M3) (t;27+T"t;4) 



IIZ 

1 -^ J ^— 3 



1 -^ 






1 -*- 



y^: 



1 -• ] ■ — 3 



1 "'-^ j-3 

2 r^ J^r^^^.'^j\ ^ 



^1 






'S'4,2 
'S'4,3 



U+U+U+U+ iki - kfy 

Aj-| AjO '"-'Q '"-'4 ^ ^ ^^ ^ 

2A (t;2T"7+Mi) (t;47+T''M3) 

ktk+k+kt W^H? 

A {uiT^^ii-i+ilT^-^u^) 

k+k+kikt ikt - kf) 



ktktk^kt iK + kt) 

k+k+k+k+ {kt - kty 

A (txiT"V37+/4T"^t^2) 

\jktktktkt (kt-kt) 

i-^^3 (A:3+ - A:4+) A (txiT" 7+ 1;2) (632C"e4) 

.^t^4 6-2 ,/^t^tktkt (kt + kty 

^'—Y^--^3 {kt + kt){kt + kt)A (e*C"e3) (e^C^) 

£ -^..^^jA„-j^^ 4 ' /7+7+7+7+ (kt — kty 

Aji Ajo '^s ''-'4 

2A;3+A;4+A (e^C^e^) (e3C"^e4) 



5' 



7,2 



^ktktk+kt (kt + A;2+)^ 



Y 1^1 1^2 "^3 "^4 

y "-1 1^2 "'S "^4 
-^^ ^^ '^ Jktk^k^kf 



Table 3: The seagull interaction in terms of Dirac spinors. The matrix elements 
Sn,j are displayed on the right, the corresponding (energy) graphs on the left. All 
matrix elements are proportional to A = g'^5{ki +k2\k^ + k4)5'''^^ {k±^i + k±^2 \ ^±,3 + 



k±,4), with g = g/J2{2'K)^. For periodic boundary conditions one uses g = g/yVL. 
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and P are second-quantized Heisenberg operators. Unfortunately, this method is 
comphcated by its non-covariant structure and the necessity to first understand its 
comphcated vacuum eigenstate over all space and time. The presence of the square 
root operator presents severe mathematical difficulties. Even if these problems 
could be solved, the eigensolution is only determined in its rest system (P = 0); 
determining the boosted wave function is as complicated as diagonalizing H itself. 
This is why instant form wave function cannot be applied in practice to scattering 
problems. Structure functions for example cannot be calculated. 

In principle, the front form approach works in the same way. One aims at 
solving the Hamiltonian eigenvalue problem 

^l^)= p+^ l^). (14) 

which for several reasons is easier: Contrary to Pz the operator P+ is positive, 
having only positive eigenvalues. The square-root operator is absent. The boost 
operators are kinematic. Having determined the wave function in a particular 
frame with fixed total momenta P+ and P± the kinematic boost operators allow 
to covariantly transcribe to any other frame. In fact, as discussed below, one can 
formulate the theory frame-independently. 

The ket |\1/) can be calculated in terms of a complete set of functions \fi) or 

d[fj] l/i) (^1 = Y^ d[fin] |/i„) (/in I = 1- 
•' n •' 

The transformation between the complete set of eigenstates |\1/) and the com- 
plete set of basis states |/i„) are then (/i„|\l/) and usually called the wavefunctions 
\l/„//i)(/i) = (/i„|\l/). In addition to the quantum numbers of the Lorentz group, the 
eigenfunction is labeled by quantum numbers like charge, parity, or baryon number 
which specify a particular hadron h, thus 



l"^) = Y J d[fln] |/i„)^„//,(/i). 



One constructs the complete basis of Fock states |/i„) in the usual way by applying 
products of free field creation operators to the vacuum state |0): 



n = 

n = 1 

n = 2 

n = 3 



0) , 

qq■.kt,h^,\^) = b^iqi)d\q2) |0) , 

qqg : kt,k_Li,Xi) = h\qi) d\q2) a) {q-^) |0), 

gg ■.kt,k±i,Xi'^ = a^{qi)a^{q2) |0), 

; ; |o) . 



The operators b^{q), d\q) and a){q) create bare leptons (electrons or quarks), 
bare anti-leptons (positrons or antiquarks) and bare vector bosons (photons or 
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gluons). All of these particles are 'on-shelF, {k^k^)i = mf. The various Fock- 
space classes are conveniently labeled with a running index n. Each Fock state 
lUn) = \n : k^, k±i, Aj) is an eigenstate of P"*" and P± and the free part of the energy 
Pq , with eigenvalues 

i(^n i^n idn i 

The free invariant mass square of a Fock-state is Mq = {pi + P2 + ■ ■ ■ +PnJ^, thus 

M'o = Po'Po,, = P^Po -Pl = P^ (e ^^^^) - Pi- (15) 

The Fock and the physical vacuum have eigenvalues 0. 

The restriction to A;^ > is a key difference between light-cone quantization and 
ordinary equal-time quantization. In equal-time quantization, the state of a parton 
is specified by its ordinary three-momentum k = {kx, ky, kz). Since each component 
of k can be either positive or negative, there exist zero total momentum Fock states 
of arbitrary particle number, and these will mix with the zero-particle state to build 
up the ground state, the physical vacuum. However, in light-cone quantization each 
of the particles forming a zero-momentum state must have vanishingly small fc+. 
The free or Fock space vacuum |0) is then an exact eigenstate of the full front form 
Hamiltonian H, in stark contrast to the quantization at equal usual-time. However, 
the vacuum in QCD is undoubtedly more complicated due to the possibility of 
color-singlet states with P"*" = built on zero-mode massless gluon quanta, but 
the physical vacuum in the front form is still far simpler than in the usual instant 
form. 

Since kf > and P^ > 0, one can define boost-invariant longitudinal momen- 
tum fractions 

k^ 
Xi = -^ , with < a;j < 1, 

and boost-invariant intrinsic transverse momenta k_j_i Their values are constrained, 

^Xi = l and ^k±i = 0, (16) 

corresponding to the intrinsic frame P± = 0. All particles in a Fock state |/in) = 
\n : Xi, k±i, Aj) have a boosted four-momentum 

Pi = [P ,P±,P )i = XiP^, k±i + XiP±, — 

x,P+ 



The free invariant mass square of the Fock state, Eq.(15), is therefore 



j^2 ^y ( m^i + {k±i + XiPi)A _ ^2 _ f m^ + kl 

i<^n \ -^i I i6n \ -^ 
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as a direct consequence of the transverse boost properties. 

The phase-space differential d[iin] depends on how one normahzes the single 
particle states. In the convention where commutators are normalized to a Dirac- 
delta function, the phase space integration is 

/ d[nn] ■■■ = X! / [dxi(fk±i'^ . . . , with 

dxicfk±i = 5(1 — X! ^i)^ ( X! ^J-i) '^^1 • • • '^^Nn d'^k±i . . . d'^k±N„, 

where Nn is the number of particles in Fock state /i„. The additional Dirac 6- 
functions account for the constraints (p^G]). The eigenvalue equation (|14D stands 
then for an infinite set of coupled integral equations 

^ jWn'] {n ■■ Xi, kj_i, Xi\H\n : x-, k^i, A-) ^„7h(a;-, k^i, A-) 

n' 

= p^ ^ '^n/h{xi,k±,Xi), for n = l,...,oo. (17) 

Since P+ and P± are diagonal operators one can rewrite this equation as 

Y^JWJ {n : Xi, h„ \\HP+ - P^\n' : x[, k^^, K) ^n'/d<, k^,, K) 

n' 
= M^^r./h{xAA^)■ (18) 

It is therefore possible to define a 'light-cone Hamiltonian' as the operator 

Hlc = HP+ - Pi = P^P^, (19) 

so that its eigenvalues correspond to the invariant mass spectrum Mi of the theory. 
Eq.(p!8D thus stands for 

Hlc\'^) = M^^) , (20) 

in analogy to Eq. ([T^) . 

The Lorentz invariance of H^c and the boost invariance of the wave functions 
reflects the fact that the boost operators are kinematical. In fact one can boost 
the system to an 'intrinsic frame' in which the transversal momentum vanishes 
Pl = 0, thus Hlc = P^P^- The transformation to an arbitrary frame with finite 
values of Pl is then trivially performed. Consider a pion in QCD with momentum 
P = {P~^, P±) as an example. It is described by 

00 „ 
: P) = Y^ d[fin] n : XiP^ , k±i + XiP±, A,) ^n/7r(a;j, k±i, Xi), 



n=l' 



where the sum is over all Fock space sectors of Eq.(|T^). The ability to specify 
wavefunctions simultaneously in any frame is a special feature of light-cone quanti- 
zation. The light-cone wavefunctions "^n/n do not depend on the total momentum, 
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Figure 2: The Haniiltonian matrix for a meson. The matrix elements are repre- 
sented by energy diagrams. Within each block they are all of the same type: either 
vertex, fork or seagull diagrams. Zero matrices are denoted by a dot (■). The single 
gluon is absent since it cannot be color neutral. 



since Xi is the longitudinal momentum fraction carried by the i*^ parton and k±i is 
its momentum "transverse" to the direction of the meson; both of these are frame- 
independent quantities. They are the probability amplitudes to find a Fock state 
of bare particles in the physical pion. But given these light-cone wavefunctions 
'^n/h{xi,k±i, Xi), one can compute any hadronic quantity by convolution with the 
appropriate quark and gluon matrix elements, see for example 0. 

In addressing to solve Eg . (|20|) one faces several major difficulties, among them 
that the above equations are ill-defined for very large values of the transversal mo- 
menta ('ultraviolet singularities') and for values of the longitudinal momenta close 
to the endpoints x ~ or x ~ 1 ('endpoint singularities'). One has to introduce 
cut-offs A to regulate the theory in some convenient way. Subsequently one has 
to remove the cut-off dependence by renormalization group analysis. Renormali- 
zation theory is known however only for perturbation theory, for example when 
calculating Feynman scattering amplitudes in a certain order. Renormalization 
theory is not available for the bound-state problem. - But even if one has found a 
convenient regularization scheme, one faces the large (infinite) number of coupled 
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integral equations. Their nature resides in the comphcated structure of the kernel 

{n\H\n) = {n : Xi, k±i, Xi\H\n : x'j, k^^, A^). 

In analyzing the structure of this very complicated many-body problem, as done 
in Figure ^, one realizes that most of its matrix elements are zero by nature of 
the operator structure of the Hamiltonian. The Hamiltonian is zero for all sectors 
whose particle number difference is larger than 2. As an example, consider the 
fork interaction given in Eq.(|TBp particularly its matrix element F^. It scatters a 
quark into an other momentum state and distroys two gluons. In the block matrix 
element {1\H\6) = {qq\F^\qqgg) one has many ways to realize that, without that 
the anti-quark q changes its momentum. In Fig. ^ all matrix elements in the block 
(l|if|6) are therefore represented diagrammatically by the same energy diagram 
as in Table |^. Usually, the typical seagull matrix elements occur for the diagonal 
blocs of this matrix, with one exception: In the graph 5*6 of Table ^ a gg-pair is 
annihilated and scatteredg instantaneously into a (^gf-pair. Correspondingly, not all 
entries of the block matrix element (5|if |5) = {qqg\Se\gg g) can vanish. In order 
to cope with the formidable many-body problem imposed by Eq. (|2D|) , the method 
of discretized light-cone quantization is very useful. 

3.1 Perturbation theory in the front form 

Let us devote a section to the perturbative treatment of front-form gauge theory. 
Light-cone perturbation theory is really Hamiltonian perturbation theory, and we 
give the complete set of rules which are the analogues of the Feynman rules. We 
shall demonstrate in a selected example, that one gets the same covariant and 
gauge-invariant scattering amplitude as in Feynman theory, see also 0]. 

The Green's functions Gfi{x^) are the probability amplitudes that a state start- 
ing in Fock state \i) ends up in Fock state |/) at a later time x~^ 

(/ I G(x+) I t) = (/|e-^^+-"K) = z|^e--"(/|G(e)|z). 
Its Fourier transform (/ | G{e) \ i) is called the resolvent of the Hamiltonian if, i.e. 



e-Ho-U + tO^ 



Separating the Hamiltonian H = Hq + U into a free part Hq and an interaction U, 
one can expand the resolvent into the series 

(/|G(e)K) = (/|G(e) + G(e)f/G(e) + G{e)UG{e)UG{e) + ... \t). 

The rules for x'^-ordered perturbation theory follow when the resolvent of the free 
Hamiltonian G{e) = l/(e — Hq + iO+) is replaced by its spectral decomposition 

G{e) = Y^ Gn{e), G„(e) = / (i[/i„] n : kf, k±i, Aj^ — (n : kf, k±i, Aj , (21) 



n=l 
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with the energy denominator A„ = e — J2i£n ((^± +?7i^)/A;+ j +iO^. The sum runs 
over all Fock states n intermediate between two interactions U. 

To calculate then {f\{e)\i) perturbatively, all x"'"-ordered diagrams must be con- 
sidered, the contribution from each graph computed according to the rules of old- 
fashioned Hamiltonian perturbation theory [jl3|, |T4| : 



1. Draw all topologically distinct ^'''-ordered diagrams. 

2. Assign to each line a single particle momentum k^, a helicity A, as well as color 
and flavor. With fermions (electrons or quark) associate a spinor Ua{k,X), 
with antifermions Va{k,X), and with vector bosons (photons or gluons) a 
polarization vector e^lk, A). 

3. For each vertex include the matrix element {n\V\n') between Fock state n 
and n' as given in Table |l[ 

4. For each intermediate state there is an energy denominator 1/A„ in which 
e = PgTin is the incident free light-cone energy. 

5. To account for three- momentum conservation include for each intermediate 
state the delta-functions 6(P^ — Y.i kf\ and 5''^^ \P±_ — Y,i k±i]. 

6. Sum over all internal helicities (and colors for gauge theories) and integrate 
over each internal k with the weight J (Pk±dk^9{k^){2n)^^^'^. 

7. Include a factor —1 for each closed fermion loop, for each fermion line that 
both begins and ends in the initial state, and for each diagram in which 
fermion lines are interchanged in either of the initial or final states. 

8. Imagine that every internal line is a sum of a 'dynamic' and an 'instantaneous' 
line, and draw all diagrams with 1,2,3,... instantaneous lines. 

9. Two consecutive instantaneous interactions give a vanishing contribution. 

10. For each instantaneous line insert a factor (n|iy|n'), with the matrix element 
given in Table ^ and ||. 

The light-cone Fock state representation can thus be used advantageously in per- 
turbation theory. The sum over intermediate Fock states is equivalent to summing 
all x"*"— ordered diagrams and integrating over the transverse momentum and light- 
cone fractions x. Because of the restriction to positive x, diagrams corresponding 
to vacuum fluctuations or those containing backward-moving lines are eliminated. 
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3.2 The gg-scattering amplitude 

The simplest application of the above rules is the calculation of the electron-niuon 
scattering amplitude to lowest non-trivial order. But the quark-antiquark scatter- 
ing is only marginally more difficult. We thus imagine an initial (g, g)-pair with 
different flavors / 7^ / to be scattered off each other by exchanging a gluon. 

Let us treat this problem as a pedagogical example to demonstrate the rules. 
Rule 1: There are two time-ordered diagrams associated with this process. In the 
flrst one the gluon is emitted by the quark and absorbed by the antiquark, and 
in the second it is emitted by the antiquark and absorbed by the quark. For the 
flrst diagram, we assign the momenta required in rule 2 by giving explicitly the 
initial Fock state |g, q) = -^ Ylc=i KfiK^ \)dcf{^qi \)\^)- Note that it is invariant 
under SU{n^). The flnal Fock state is \q',q') = ^ Ec=i &I/(fcg, A;)4/(fcg-, A;-)|0). 
The intermediate state 



W,(l,9) 




-■- c=l c'=l a=l 



has 'a gluon in flight'. Since the gluons longitudinal momentum is positive, the 
diagram allows only for k'^ < k^. Rule 3 requires at each vertex the factors 

{QyQ]^ IQ ,(1,9/ = 3"W^; ; ; ; -, US) 

(g',g,f7|\/|g',g') = -f-j\H=. ^ ; ^ == -, (24) 




2ktpkj^2k',^ 

respectively, wich are found by means of Table |I[ Working with color neutral 
Fock states, all color structure reduces to an overall factor C, with C^ = (n^ — 
l)/2nc. Rule 4 requires the energy denominator I/A3. It is useful to work the 
four-momentum transfers of the quark Q^ = —{kg — k'^)"^ = kg {kg + k'^ — kg)~ . 

9 

The anti-quark has Q = {kg — k'^y = k^{kg + kg — k'^)~ . With the initial energy 
e = P^ = {kg + kg)+ = ^{kg + kq)~ , the energy denominator becomes then 

A3 = {kg + kg)- - {kg + k'g + kg)' = "^ 

Rule 5 requires two Dirac-delta functions, one at each vertex, to account for conser- 
vation of three-momentum. One of them is removed by the requirement of rule 6, 
namely to integrate over all intermediate internal momenta and the other remains 
in the flnal equation (^61) . The gluon momentum is thus flxed by the external legs 
of the graph. The polarization sum over the gluon helicity gives 

, / J \ v~^ /I \ \ */] \ \ 9tfJ'Vi' ~r rigiiTj^ 

dfiuyl^g) = 2^ ^fiyl^g, ^g) ^ui'^gy ^g) = ~9fiu H 7— • 
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The null vector has the components r]^ = ij]'^ ,f]^,rj~) = (0,0i,2) and thus the 
properties rj'^ = rj^rj^ = and krj = k^ . As shown explicitly in [^ one gets for the 
second order diagram {q-,(i\VG^V\q\q') after some non-trivial steps 



{q,q\VG^V\q',q') 



2r^2 



(2vr) 



m(A;<?, Ag)7^u(A;^, A'J i \u{kq,Xq)-f^u{k'q,X'q) 



Ak^k'+ 



Q' 



4fc,7 K^ 



g^c^ \u{kg, Xqh+uiK^ K) 



q '"q 
U{kq, Xq)-f+u{k' X'g] 



{2n) 



Ak+k'+ 



{K 



Akq kg 



(25) 



The delta-functions and a step function 0(/cl^ < k'^) are omitted for simplicity. 
One proceeds with rule 8, by including the instantaneous lines. Table gives 



{q,q\S\q',q') 



g^C' [u{k,X)j+u{k',X% 
(27r)3 



[u{k,X)j+u{k',X% 



Ak+k'g^ 



iK - K' 



Akq kg 



The gg-scattering amplitude, the sum (g, q\ S + VG^V \q', q') has then the correct 
gauge-invariant result known from Feynman theory up to second order 



{q,q\S + VG^V\q\<i) 



iK-K)\/k+ktk'+k',+ 



S{P+-P'+)6^^\P^-Pl) 



q '"g "^q '"q 



(27r)3 



[uik, X)Yuik', X% [u{k, Xh,u{k', X% . (26) 



The instantaneous diagram is cancelled exactly. 



4 Discretized Light-Cone Quantization 

The infinitely many coupled integral equations in the preceeding section a very 
difficult to cope with in practice. Because of the many integrations and summa- 
tions it is very difficult to even write them down. But field theory becomes much 
more transparent when one works with periodic boundary conditions. The coupled 
integral equations become then coupled matrix equations. Since rows and colums 
of a matrix can be denumerated, one can keep track of the necessary manipula- 
tions much easier. In fact, working with periodic boundary conditions, or with 
'Discretized Light-Cone Quantization' (DLCQ), one has obtained the first total so- 
lutions to non-trivial quantum field theories in 1-1-1 dimensions. In 3-1-1 dimensions 
the method has the ambitious goal to calculate the spectra and wavefunctions of 
physical hadrons from a covariant gauge field theory. The ingredients of the method 
shall be reviewed in short in this section. 
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Figure 3: Non-relativistic many-body theory. 



4.1 The non-relativistic A-body problem in one dimension 

Let us first briefiy review the difficulties for a conventional non-relativistic many- 
body theory. One starts with a many-body Hamiltonian H = T + U. The kinetic 
energy T is usually a one-body operator and thus simple. The potential energy 
U is at least a two-body operator and thus complicated. One has solved the 
problem if one has found the eigenvalues and eigenfunctions of the Hamiltonian 
equation, H'^ = E"^. One always can expand the eigenstates in terms of products 
of single particle states {x\m) belonging to a complete set of ortho-normal functions. 
When antisymmetrized, one refers to them as 'Slater-determinants'. All Slater- 
determinants with a fixed particle number form a complete set. 

One can proceed as follows. In the first step one chooses a complete set of 
single particle wave functions. These single particle wave functions are solutions of 
an arbitrary single particle Hamiltonian. In a second step, one selects one Slater 
determinat as a reference state. All Slater determinants can be classified relative 
to this reference state as 1-particle-l-hole (1-ph) states, 2-particle-2-hole (2-ph) 
states, and so on. The Hilbert space is truncated at some level. In a third step, 
one calculates the Hamiltonian matrix within this Hilbert space. 

In Figure ^, the Hamiltonian matrix for a two-body interaction is displayed 
schematically. Most of the matrix-elements vanish, since a 2-body Hamiltonian 
changes the state by up to 2 particles. Therefore the structure of the Hamiltonian 
is a finite penta-diagonal block matrix. The dimension within a block is made finite 
by an artificial cut-off on the kinetic energy, i.e. on the single particle quantum 
numbers m. A ffnite matrix can be diagonalized on a computer. At the end one 
must verify that the physical results are reasonably insensitive to the cut-off(s) and 
other formal parameters. 

This procedure was actually carried out in one space dimension ||15[ with two 
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different sets of single-particle functions, 

(a;|m) = A^„ifm(-) exp|--(-) I and {x\m) = N^exp lim-ni . (27) 

The two sets are the eigenf unctions of the harmonic oscillator {L = h/muj) with its 
Hermite polynomials H^, and the eigenfunctions of the momentum of a free particle 
with periodic boundary conditions. Both are suitably normalized (Nm), and both 
depend parametrically on a characteristic length parameter L. The calculations 
are particularly easy for particle number 2, and for a harmonic two-body inter- 
action. The results are displayed in Figure ^, and surprisingly different. For the 
plane waves, the results converge rapidly to the exact eigenvalues E = ^,^,y, ■ ■ ■, 
as shown in the right part of the figure. Opposed to this, the results with the 
oscillator states converge extremely slowly. Obviously, the larger part of the Slater 
determinants is wasted on building up the plane wave states of center of mass mo- 
tion from the Slater determinants of oscillator wave functions. It is obvious, that 
the plane waves are superior, since they account for the symmetry of the problem, 
namely Galilean covariance. The approach was successfully applied for getting the 
exact eigenvalues and eigenfunctions for up to 30 particles. 

From these calculations, one should conclude that discretized plane waves are 
a useful tool for many-body problems, and that the generate good wavefunctions 
even for a 'confining' potential like the harmonic oscillator. 

4.2 QED in 1+1 dimension 

DLCQ had been applied first to Yukawa theory [l^ in 1-space and 1-time dimen- 
sions followed by an application to QED |T^ and to QCD ||T^, but the advantages 
of working with periodic boundary conditions in the front form, particularly when 
discussing the 'zero modes' in the 0^-theory, had been noted also by Maskawa and 



Yamawaki [^] in 1976. 

In one space dimension are no rotations — hence no spin. The Dirac spinor has 
two components and the Dirac matrices are 2 by 2 matrices. The gauge field A'^ 
field has two components. One of them is eliminated by fixing the gauge, and the 
other determined by Gauss' law. The gauge field carries no dynamical degree of 
freedom. The theory confines quarks (or electrons) because the Poisson equation 
in 1 space dimension gives rise to a linearly rising potential. 

Quantum electrodynamics in 1-1-1 dimension has played an important role in 
field theory because its massless version, the Schwinger model, is analytically 
soluable. The fundamental solution corresponds to a composite ee-state - the 
Schwinger-boson - with invariant mass tub = g/y/^- Note that the coupling con- 
stant g in 1-1-1 dimension has the dimension of a mass. 

Consider first the massive Schwinger model. The Lagrangian for the theory 
takes the same form as in Eq.(|^). Again one works in the light-cone gauge A'^ = 0, 
and uses the same projection operators A±. The Dirac equation decomposes again 
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Table 4: Fock-space sec- 
tors and block matrix struc- 
ture for QED in 1+1 di- 
mension. Diagonal blocs 
(D=T-I-S) refer to seagull, 
off-diagonal (F) to fork ma- 
trix elements. A zero block 
matrix is marked by (■). 



into a time-derivative 2i9+\l/-|- = r7i/5\l/_ + gA^'^^ and a constraint 229^\E'_ = 
m/3\l/+. One is left with only one independent field, \E'+, which is canonically 
quantized at x"*" = 0, 



{'^+{x-,x+),'^\{y~,y+)} 



x+=y+=0 



A+6{x" -y^). 



f28) 



The formalism gets somewhat simplified if one uses the chiral representation ||18 
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The A± = (1 ± 7°7^)/2 are then diagonal and project on the chiral components 
The (light-cone) momentum and energy operators become 
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i^£l^d.'¥,v'^^,„,. 



One can expand ^ji (or \E'_|_) with periodic boundary conditions ||T3, but anti- 
periodic boundary conditions ||18[ avoid the zero mode: 



^i?(a;~ 



1 



2L 



Y. i^nC 



+ cit e^¥- 



The creation and destruction operators obey 1 6jj, 6^ !• = I S^, dm \ 
tent with Eg . (p8[) . One redefines units by 



6n,m, consis- 



^^-f*- n-lp^o + ^C/ 



7T 



4 DISCRETIZED LIGHT-CONE QUANTIZATION 



28 



1 


1. 7 
2' 2 




2 


3. 5 




2 


2' 2 




4 


7. 1 
2' 2 




5 


1 3. 1 


3 


2' 2' 2' 


2 



K 


1 qq 


2qq 


3qq Aqq 


Total 


1 


1 


- 


- 


1 


4 


4 


1 


- 


5 


9 


9 


20 


1 


30 


16 


16 


140 


74 1 


231 



Table 5: The 5 Fock states for K = A. Table 6: Number of Fock states versus 

harmonic resolution. 



The length L drops out in i^LC = P~^P ■ Inserting the above fields gives 

OO OO -I 

K = Yl n{b% + did^), Ho= J2 -{blK + didr,), 



1 3 

"2 '2 ' 



1 3 
"2'2' 



n 



U 



■.U : 



CXJ T 

■.U:+ Yl ~ {^li^n + didn) , with I„ 
1 ^ ^ 



1 3 

"2'2' 



1 

'2n 



E 

m=l 



1 

2 1 



m,^ 



E 



ni,n2,n3,n4 



5{ni + n2\n3 + n^) 
2(ni - 723)2 



(6162^364 + (i|(i2C^3d 



+ 



+ 



Y '^('^i + ^2|n3 + ^4) I - — — 

Il,n2,n3,n4 \ l"! 



E 



ni,n2,n3,n4 



6{ni\n2 + ns + 714) 
{rii - n^y 



- n^Y (ni + 722)2 ^ 
\h\h2h3di + d\d2dsb4^ + h.c.j 



bldlb^d^ 



The symbols 6{n\m) = 5n,m are Kronecker delta's and 61 = 6„^. In 1+1 dimensions 
it is very important to keep the 'self-induced inertias' /„ from the normal ordering. 
They are needed to cancel the infrared singularity in the interaction term in the 
continuum limit. 

The next step is to actually solve the equations of motions in the discretized 
space. Typically one proceeds as follows: One constructs the Fock space 

I/in) = |gi,...,g„;gi,...,g„) = b\^, . . . ,b\j\^, . . . , d\j^Q) , 

in the same way as above in Eq.(|TBp, and arranges it in denumerated Fock space 

sectors, as illustrated in Table |. Each Fock state must be an eigenstate to P^ and 

thus of the harmonic resolution K [0, with eigenvalue K = J^i^/i^ '^j- One selects 

now one value for K and constructs all Fock states. Both the number of Fock 

space sectors and the number of Fock states within each sector are finite for every 

finite K, due to the positivity condition on the light-cone momenta. For K = 1 

one has only one Fock-space class with one Fock state |^; ^) = b\di\0), ioi K = A 

22 
one has the 5 Fock states given in Table ^. The numbers of Fock states increase 

with given K, as shown in Table 0, but less than exponentially due to the exclusion 
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Figure 4: Spectra and wavefunctions in 1+1 dimensions, taken from |Ij, \T8 



principle. Next, one calculates the matrix elements of H = P^ . In the last step 
one diagonalizes H. Any of its eigenvalues E{K) depends on K and corresponds 
to an invariant mass M'^{K) = P^P~ = KE[K). Notice that one gets a spectrum 
of invariant mass-squares for any value of K. 

The eigenvalue spectrum of QEDl + 1 was given first by Eller 



17| , for periodic 



boundary conditions on the fermion fields. The plot likeon the left side in Figure ^ 
was calculated p2| with anti-periodic boundary conditions. It shows the full mass 
spectrum of QED in the charge zero sector for all values of the coupling constant 
and the fermion mass, parametrized by A = (1 + n{m/gY)~2. The eigenvalues 
Mi are plotted in units where the mass of the lowest 'positronium' state has the 
numerical value 1. All states with M > 2 are unbound. The plot includes the free 
case X = {g = 0) and the the Schwinger limit M = 1 for A = 1 {m, = 0), where 
DLCQ generates the exact eigenvalue. - The lower left part of the figure illustrates 
the following point. The rich complexity of the spectrum allows for multi-particle 
Fock states at the same invariant mass as the 'simple gg-states' shown in the figure 
as the '2 particle sector'. The spectrum includes not only the simple bound state 
spectrum, but also the associated discretized continuum of the same particles in 
relative motion. One can identify the simple bound states as two quarks connected 
by a confining string as displayed in the figure. The smallest residual interaction 
mixes the simple configuration with the large number of 'continuum states' at the 
same mass. The few simple states have a much smaller statistical weight, and it 
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looks as if the long string 'breaks' into several pieces of smaller strings. Loosely 
speaking one can interpret such a process as the decay of an excited pion into 
multi-pion configurations tt* -^ tttttt. 

In 1 + 1 dimensions quantum electro dynamics |]T^ and quantum chromo dynam- 
ics |l3| show many similarities, both from the technical and from the phenomeno- 



logical point of view. In the right part of Figure ^ some of the results of Hornbostel 



T8|| on the spectrum and the wavefunctions for QCD are displayed. Fock states 
in non-abelian gauge theory SU(N) can be made color singlets for any order of 
the gauge group and thus one can calculate mass spectra for mesons and baryons 
for almost arbitrary values of N. In the upper right part of the figure the lowest 
mass eigenvalue of a meson is given for A^ = 2, 3, 4. Lattice gauge calculations are 



available only for N = 2 and for the lowest two eigenstates [0. In general the 
agreement is very good. In the left lower part of the figure the structure function 
of a baryon is plotted versus (Bj0rken-)a; for m/g = 1.6. With DLCQ it is possible 
to calculate also higher Fock space components. As an example, the figure includes 
the probability distribution to find a quark in a ggggg -state. 

4.3 Fermion condensates and the small mass limit 

Based on the low energy theorems from times prior to QCD, it is believed that 
the square of the pion mass is linear in the quark mass m for sufficiently small 
m. The proportionality constant has to have a dimension of mass, and since 
there is no other scale in the problem except the quark condensates in the vacuum 
(0|\E'\E'|0), one believes that the square of the pion mass is approximatively given 
by rri^ ~ 2(0|\I'\&|0)m, a theorem which was succesful in many phenomenological 
applications. 

The Schwinger model has played an important role as a paradigm for our un- 
derstanding of hadronic physics. Among other aspects it has the desired feature 
that the invariant mass square of the ee-boson is linear in the electron mass m 

M^ = m| + 2m(0|^*|0), with (0|^^|0) = e^ m^, (29) 

where ms = QJ \fTi is the invariant mass in the Schwinger limit. Euler's constant 
7 was understood as a signal for non-perturbative physics. In his analysis of the 



Schwinger model, Bergknoff |^ showed that the take-off from the Schwinger limit 
obeys 

(0|^^|0) = ^ vfiB. (30) 

The result of this 'chiral perturbation theory' to first order, ^ ~ 1.81, is numeri- 
cally very close to e"^ ~ 1.78. 

Its is actually quite easy to derive the Bergknoff equation from DLCQ ||21[|. For 
a fixed resolution K^ a gg-state is fixed by the quantum number n of the electron. 
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Figure 5: The lowest mass eigenvalue 
of QEDi+i in units ofms is plotted ver- 
sus the ferniion mass ruf in units of the 
coupling constant g, as calculated with 
'conventional DLCQ' (K = 16 j. Taken 
from Ref p^/ . See also discussion in the 
text. 
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Figure 6: The lowest mass eigenvalue 
of QEDi+i in units ofms is plotted ver- 
sus the fermion mass ruf in units of the 
coupling constant g, as calculated with 
'improved DLCQ' (K = 20), as pro- 
posed in Ref J2^. Taken from Ref. pBjj . 
See also discussion in the text. 



The eigenvalue equation is then given by 



[K] 



M2(n|^) = J2 in\HLc\n') {n'\^), for 



n 



1 3 

2'2' 



[Kl 



with [K] = K — ^. Using the Hamiltonian matrix elements are given above, one 
introduces x = p'^ /P^ = n/K and goes to the continuum limit. After a few steps 



2T[| one gets the integral equation of Bergknoff P0 | 

1 



M2(x|^) 



m, 



x{l — x) 



(x|\I') + m| / dx'{x'\'^) +m%-h dx 



,(x|^) - (x'l*) 



(x — x') 



.l\2 



For 771 = 0, the solution (x|\l') = 1 has the eigenvalue M^ = m?^, the Schwinger 
boson. 

Eqs.(pUD and ( pil| ) state that the mass-squared of the Schwinger boson is linear 
in the fermion mass m, in the limit when m goes to zero. For a long time, this result 
was in conflict with the explicit DLCQ-calulations [|T^, |22|, ^. The shortcoming 
was taken as a hint ||2^ that light-cone quantization was in failure because of the 
trivial vacuum structure which does not allow for condensates. 
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The most recent results by Vollinger [^ for DLCQ are displayed in Figure |. 



They are compared with the lattice calculations of Crewther and Hamer [^ and 



with chiral perturbation theory up to second order [0]. The lattice results and 
DLCQ show some discrepancy for very small m which however fades away for 
larger m. On the other hand, the lattice results and chiral perturbation theory 
agree perfectly at small and deviate for large m, where chiral perturbation theory 
is not suposed to work. The figure show aslo that a correct inclusion of the (gauge 
field) zero modes has no significant impact. The re-solution of this puzzle came 



by van de Sande |^. He realized that conventional DLCQ has difficulties to 
reproduce the wave function (x|\l/) near the endpoints x — ^ and x —>■ 1 very close 
to the Schwinger limit. His 'improved DLCQ' accounts for that, and indeed when 
properly included as shown in Figure ^ all discrepancies fade away. 

What should be learned from this exercise is that not everything what is called 
a 'vacuum condensate' in the literature deserves this name in a physical sense: In 
naive light-cone quantization particularly DLCQ the vacuum is trivial and can not 
have condensates. 

4.4 $^ in 1+1 dim's: Zero modes and phase transitions 

The naive front-form vacuum is simple. However, one commonly associates im- 
portant long range properties of a field theory with the vacuum like spontaneous 
symmetry breaking, the Goldstone pion, or color confinement. If one cannot asso- 
ciate long range phenomena with the vacuum state itself, then the only alternative 
is the zero momentum components of the field, the 'zero modes'. In some cases, 
the zero mode operator is not an independent degree of freedom but obeys a con- 
straint equation. Consequently, it is a complicated operator-valued function of all 
the other modes of the field. Zero modes of this type have been investigated first 
by Maskawa and Yamawaki as early as in 1976 ||2^. An analysis of the zero mode 



constraint equation for (1-|-1) -dimensional 0^ field theory, by van de Sande and 
Pinsky [^ , shows how spontaneous symmetry breaking occurs within the context 



of this model. 

The model represents a new paradigm for spontaneous symmetry breaking and 
shall be reviewed shortly. The Lagrangian in one space and one time dimension 
is £ = 9_|_09_(^ — ^0^ — li^^- Imposing periodic boundary conditions with a 
length parameter d = 2L gives 0(x) = 4^ SJ^-oo qn{x^)e^^"'^ ■ The field integral 
T,n = J dx~ (j){x)^ — {zero modes) is convenient to discuss the problem and becomes 

'^n = — 2^ ?ii?i2 ■■•?«« ^n+i2 + -..+Jn,0- 

Also the canonical Hamiltonian can be disentangled into E„ and the zero mode go 
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Figure 7: The vacuum expectation value 
of (f), /o = \/47r(O|0|O) is plotted versus 
g = 247r/i^/A, which is the inverse of the cou- 
pling constant A. Taken from ^2^/ . — In the 
front form, the vacuum state |0) is simple, 
but the operator (p, or oq, is complicated. In 
the conventional instant form, the vacuum 
|0) is complicated but the operator cf) is sim- 
ple. 



One can apply canonical quantization. Following the Dirac-Bergman prescription, 
described in 0, one identifies first-class constraints which define the conjugate 
momenta = Pn — ik^^Q-n ,, where [qm,Pn] = Sn,m/'^ and m, n 7^ 0. The secondary- 
constraint 







(31) 



determines the zero mode go- This result can also be obtained by integrating the 
equations of motion, 9+9_0 + /i^0 = |y0^. To quantize the system one replaces the 
Dirac bracket by a commutator. One must choose a regularization and an operator- 
ordering prescription in order to make the system well-defined. One begins by 
defining creation and annihilation operators a], and a^. 



Qk 



d 



A7T\k\ 



a-k 



a-k 



-k ■> 



k^O 



which satisfy the usual commutation relations ak,ai = S^^i- Likewise, one defines 

the zero mode operator go = aoJd/An. In the quantum case, one normal orders 
the operator S„. One redefines this the fields in terms of operators 



[X 



Vd\ ^,n^ 



^^"^" + ale'^"^" 



and notes that they obey the canonical commutation relations 



(x+,x ),d-(p{y^,y ) 



6{x -y ), 



(32) 



(33) 



for a boson field in the front form, see also Eq. (p8|) . 

The solution of the constraint Eq. (|3lD is very difficult. Van de Sande and Pinsky 



| 27[] have shown that the zero mode a^ could acquire finite values, 

(O|0|O) = ao ^ 0, 
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depending on the coupling constant. 

One finds the following general behavior: for small coupling (large g, where 
g oc 1/coupling) the constraint equation has a single solution and the field has no 
vacuum expectation value (VEV). As one increases the coupling (decreases g) to the 
"critical coupling" (^critical, two additional solutions which give the field a nonzero 
VEV appear. These solutions differ only infinitesimally from the first solution near 
the critical coupling, indicating the presence of a second order phase transition. 
Above the critical coupling {g < (^critical), there are three solutions: one with zero 
VEV, the "unbroken phase," and two with nonzero VEV, the "broken phase". The 
"critical curves" shown in Figure ^, is a plot of the VEV as a function of g. 

5 DLCQ in 3+1 dimensions 

Periodic boundary conditions on C can be realized by periodic boundary conditions 
on the vector potentials A^ and anti-periodic boundary conditions on the spinor 
fields, since C is bilinear in the \&q. In momentum representation one expands 
these fields into plane wave states e"*^*"^^, and satisfies the boundary conditions by 
discretized momenta 

{ju, with n = i, |, . . . , oo for fermions, 
jTi, with n = 1, 2, . . . , oo for bosons, 

vr ^ 
and p± = -— nj_, with n^;, n^^ = 0, ±1, ±2, . . . , ±cxd for both. 

As an expense, one has to introduce two artificial length parameters, L and L±. 
They also define the normalization volume Q = 2L{2Lj_y. More explicitly, the free 
fields are expanded as 

V * ' q \ P 

and A,{x) = ^Y.^ (se/.(P, A)e-P^' + a\el{p, \)d^^) , (34) 

particularly for the two transverse vector potentials A^ = A]_, [i = 1,2). As 
above, each single particle state 'g' is specified by six quantum numbers, the three 
discrete momenta n, n^jUy, helicity, color and flavor. The creation and destruction 
operators like aj and a^ create and destroy single particle states q, and obey (anti-) 
commutation relations like 

[ag,al,] = {bg,bl,} = {dg,dl,} = 6g,q>. 

The Kronecker symbol is unity only if all six quantum numbers coincide. 

Inserting the free fields of Eq.(|3^) into the Hamiltonian, one performs the 
space-like integrations and ends up with the light-cone energy-momenta P^ = 
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P'^{aq,a}j,bq,b]j,dq,d]j) as operators acting in Fock space. The spatial components 
of P^ are simple and diagonal, the temporal component P~ is complicated and 
off-diagonal. The integrals over the coordinates x~ are conveiently expressed in 
Kronecker delta functions 

+L +L 

-L -L 

and correspondingly for the transversal integrations. In the tables given above they 
appear typically in the overall factor 

A(gi;g2,g3,g4) = ^ S{kt\k+ + k+ + kt) S^^\k^i\k^2 + A;±3 + k±4). 

5.1 Retrieving the continuum formulation 

The continuum formulation of the Hamiltonian problem in gauge field theory with 
its endless multiple integrals is usually cumbersome and untransparent. In DLCQ, 
the continuum limit corresponds to harmonic resolution K ^ oo. The compactified 
formulation with its simple multiple sums is straightforward. The key relation is 
the connection between sums and integrals 



Jdk+fik-', k^)^^^Y. /(^^' k^ 

Jd% fik+, 4) ^^ |J- E /(^^' ^j- 



Combined they yield 

fdkU'k^ f{k^, h) ^^ ^^ Y. /(^^' k^)- 

Similarly, Dirac delta and Kronecker delta functions are related by 

6{k^) 6^'\k) ^^ ^^ S{k^\0) 6^'\k\0). 

Because of that, in order to satisfy the respective commutation relations, one must 
modify also the creation and destruction operators. Denoting the single boson 
operators in the continuum by a and in the discretized case by a, they must be 
related by 



«(?) ^^ \|7U7^Z^^ «r 



and correspondingly for fermion operators. Of course, one has formally to replace 
sums by integrals, Kronecker delta by Dirac delta functions, and single particle 
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operators by their tilded versions. In practice, it suffices to replace the tilded 
couphng constant 

-2 ^^ u ~2 9^ 



r = ^ by r 



2n ' ^ 4(27r)3 

in order to convert the discretized expressions in Tables |I|-^ to the continuum 
formulation. 

The DLCQ method can be considered a general framework for solving problems 
such as relativistic many-body theories or approximate models. The general pro- 
cedure is: (1) Phrase the physics problem in DLCQ; (2) Apply approximation and 
simplifications; (3) Derive the final result; (4) At the end convert th so obtained 
expressions to the continuum. 

5.2 Fock-space and vertex regular izat ion 

The finite number of Fock-space sectors is a consequence of the positivity of the 
longitudinal light-cone momentum p"*". The transversal momenta px can take either 
sign, and the number of Fock states within each sector can be arbitrarily large. In 
order to face a finite dimensional Hamiltonian matrix one must have a finite number 
of Fock states, and this is achieved by Fock space regularization: Following Lepage 
and Brodsky [jl4|, a Fock state with n particles is included only if its free invariant 



mass does not exceed a certain threshold 



{Pi+P2 + ■■■ Pnf -{mi+m2 + ... m„f < A^. 

The sum extends over all n particles in a Fock state. The lowest possible value of 
Mq is taken when all particles are at rest relative to each other, i.e. {Mq)^.^ = 
(mi -|- 777.2 + • • • rrin) . This frozen invariant mass should be removed from the cut- 
off The mass scale Aq is a Lorentz scalar and one of the parameters of the theory. 
However, it was not realized in the past f^, that Fock-space regularization is 
almost irrelevant in the continuum theory. Vertex regularization seem to be a better 
alternative. At each vertex, a particle with four-momentum p^ is scattered into two 
particles with respective four-momentum pi and P2. In order to avoid potential 
singularities one can regulate the interaction by setting the matrix element to zero 
if the off-shell mass {pi +^2)^ exceeds a certain scale A. The condition 

R{A) = ({pi + P2)' - (mi + 7772)2 - A^) (35) 



will be referred to as the sharp cut-off for vertex regularization. 

The dependence on the regularization parameter A must be removed by the 
renormalization group. Renormalization looks like a terrible problem in the context 
of non-perturbative theory, but with the above regularization scheme it could be 
simple in principle: The eigenvalues may not depend on the regulator scale(s) A. 
To require this is easier than to find a practical realization. 
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6 The many-body problem in gauge theory 

In principle one proceeds in 3+1 like in 1+1 dimensions: One selects a particular 
value of the harmonic resolution K and the cut-off A and diagonalizes the finite 
dimensional Hamiltonian matrix by numerical methods. But the bottle neck of any 
Hamiltonian approach is that the dimension of the Hamiltonian matrix increases 
exponentially fast with the cut-off. As a concrete example consider the matrix 
structure as given ioi K = 4 in Fig. ^ Suppose the regularization procedure 
allows for 10 discrete momentum states in each direction. For every single particle 
one has about 10^ possibilities to define a momentum state. A Fock-space sector 
with n particles and fixed total momentum has then roughly 10""^ different Fock 
states. Sector 13 in Fig. ^ alone, with its 8 particles, has thus about 10^^ Fock 
states. Chemists are able to handle matrices with some 10^ dimensions, but 10^^ 
dimensions exceeds the calculational capacity of any computer in the foreseeable 
future. 

For 3+1 dimensions one is thus confronted with a a similar problem as in 
conventional many-body physics, displayed in Fig. ^. One has to diagonalize finite 
matrices with exponentially large dimensions (typically > 10^). In fact, the problem 
in quantum field theory is even more difficult since the particle number is unlimited. 
One needs an effective interaction which acts in smaller matrix spaces and which 
has a well defined relation to the full interaction. One needs it also for the physical 
understanding. The effective interaction between two electrons, for example, is the 
Coulomb interaction, at least to lowest order of approximation. 

The goal is therefore to develop an exact effective interaction between a quark 
and an anti-quark in a meson. To achieve this goal, the Hamiltonian DLCQ-matrix 
is discussed in terms of block matrices, since the Fock-space sectors appear quite 
naturally in a gauge theory, see also Eg . (pIsD . Each Fock-space sector has a finite 
number of Fock states, which are kept track of collectively. Eq.([20|) is therefore 
rewritten as a block matrix equation, with H = H^c and E = M^, 

N 

Y, {i\H\j) {il'il) = E {i]"^) for alH = l,2,...,iV. (36) 

i=i 

Rows and columns are denumerated in the same convention as in Table |^. As to 
be shown, it can be mapped identically on a matrix equation which acts only in 
sector |1). Once this is achieved, one can go to the continuum limit. 

6.1 The approach of Tamm and Dancoff 

Effective interactions are a well known tool in many-body physics 0. In field theory 
the method is known as the Tamm-Dancoff-approach, applied first by Tamm [^ 
and by Dancoff |2^, which shall be reviewed it in short. 



The rows and columns of a matrix are split into the P- and the Q-space. In 
terms of the sector numbers of Eq.(P^, P = Yl]=i \j){j\ and Q = Y.f=n+i \J){J\^ 
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where 1 < n < A^. Eq.(|3q) can thus be written as a 2 by 2 block matrix equation 

{P\H\P) {P\^) + {P\H\Q) {Q\^) = E{P\^), (37) 

{Q\H\P) {P\^) + {Q\H\Q) (Ql^i^) = E{Q\<i^). (38) 

If one can invert the quadratic matrix {Q\E — H\Q) one could express the Q-space 
in terms of the P-space wavefunction. But here is a problem: The eigenvalue E is 
unknown at this point. But one can replace it by another number, the starting point 
energy uj, which is at first a free parameter. The matrix inverse to {Q\uj — H\Q) is 
called the resolvent of the Hamiltonian in the Q-space and is denoted by Gq{uj). 
The Q-space wave function becomes then 

(Q|M/(^)) = Gq{u){Q\H\P) (P|M/), Gq{u:) = ^g^^^^^gy (39) 



Substituting it in Eq.([37|) produces an eigenvalue equation in the P-space, 

H,s{u)\P) {P\^k{uj)) = Ek{uj) |^fc(a;)), (40) 

and defines the effective P-space Hamiltonian 

H,s{uj) = H + H\Q)Gq{lo){Q\H. 

Varying u one generates a set of energy functions Ek{uj). Every solution of the 
fix-point equation 

Ek{uj) = u, 

generates one of the eigenvalues H, in fact, it generates all of them. 

If one identifies the P- with the gg-space one seems to have found the effective 
interaction which acts in the Fock space of a single quark and a single anti-quark. It 
looks as if one has mapped a difficult problem, the diagonalization of a big matrix 
onto a simpler problem, the diagonalization of a small matrix. But the price to 
pay is to invert a matrix. Matrix inversion takes about the same numerical effort 
as its diagonalization. In view of having to vary u, the numerical work is therefore 
rather larger than smaller as compared to a direct diagonalization. The advantage 
of working with a resolvent is of analytical nature to the extent that resolvents can 
be approximated systematically. The two resolvents 

defined once with and once without the off-diagonal interaction U, are identically 
related by Gq = G -\- GUGq, or by the infinite series of (Tamm-Dancoff) pertur- 
bation theory 

Gq = G + GUG + GUGUG + ..., (42) 
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Table 7: The Fock-space sectors and the Haniiltonian block matrix structure 
for QCD. Diagonal blocs are marked by D. OS-diagonal blocks are labeled by 
V , F and Sq, corresponding to vertex, fork and seagull interactions, respectively 
Zero-matrices are denoted by dots. Taken from P^ /. See also Fig. ^. 



see also Eq.(^). The free resolvent G can be obtained trivially since the kinetic 
energy T is diagonal. Conceptually, it is the same object as the one in Eg. (pT]) , 
except for two aspects: The present G acts only in the Q-space, and the starting 
point energy cj is a constant (one of the eigenvalues); the free energy e in ([21|), how- 
ever, is a function of the incoming momenta. The starting point energy not being 
a kinetic energy creates problems allover the place, since non-integrable singulari- 
ties appear in every order of (Tamm-Dancoff ) perturbation theory. Essentially two 
conclusions are possible: Either gauge theory has no bound-state solution, or the 
series in Eg. (^21) has to be resumed to all orders of perturbation theory before the 
singularities begin cancel eachother. 



In practice, Tamm and Dancoff |2^, ^ have restricted themselves to the first 
non-trivial order. In order to make things work, they have replaced the energy 
denominator with the eigenvalue uj by the energy denominator with the function e, 
as it appears in the perturbative scattering amplitudes. The same trick was applied 
in the later work with the front form |]5D|, ^ . 

Perturbation theory within a bound state problem is known to be a very difficult 
guestion. It has motivated the formal work in the next few sections. What we are 
after is the impossible, some kind of non-perturbative perturbation theory! 



6 THE MANY-BODY PROBLEM IN GAUGE THEORY 40 

6.2 The method of iterated resolvents 

The Tamm-Dancoff approach can be interpreted as the reduction of a block matrix 
dimension from 2 to 1. But having a matrix with block matrix dimension A^ = 13 
as in Table ^ it could be interpreted as the reduction from A^ -^ A^ — 1, simply by 
choosing the Q-space appropriately. But then the procedure can be iterated, one 
can reduce the block matrix dimension from A^ — 1 ^ A^ — 2, and so on until one 
arrives at 2 ^ 1. This method of 'iterated of resolvents' has certain advantages, 
which will be discussed as the formalism develops. 

First, one needs a reasonable and compact notation. One of them is to denu- 
merate the Fock space sectors as in Table |^. The full Hamiltonian will be denoted 
by if = Hn, since by the definition in Eq.(PBD it has A^ blocks. Suppose, during 
the reduction one has arrived at block matrix dimension n, with 1 < n < A^. The 
eigenvalue problem corresponding to Eq.(^0|) reads then 



, ...,n. 



Y,{t\H4uj)\j){j\^{uj)) = E{u) m{uj)), for z = 1,2 

In analogy to Eg. (|39|) , define the n-space resolvent, and get 

{n\^iu;)) = G^{uj)j:{n\H^iu;)\j){j\^iu)), G^u) = — ^ (43) 

p{ {n\uj - Hn{uj)\n} 

The effective interaction in the {n — l)-space becomes then 

Hr,.i{uj) = Hn{uj) + Hn{uj)Gn{uj)Hn{uj) (44) 

for every block matrix element. It is unpleasant but unavoidable, that the symbol 
n denotes both the block matrix dimension and the number of the last sector. Else 
one proceeds like for Tamm-Dancoff, including the fixed point equation E{uj) = uj. 
But one has achieved much more: Eq.(Hl) is a recursion relation] 



A few comments seem to be in order. The method of iterated resolvents [32, 33 
is particularly suited for gauge theory with its many zero block matrices, see Ta- 
ble §. The zero matrices remove many of the multiplications in Eq. (|43D . The 
Tamm-Dancoff procedure cannot make use of them. - Both the Tamm-Dancoff 
procedure and the iterated resolvents can be put on a computer in model studies. 
The iterated resolvents require to invert several smaller matrices insted of one big 
one. Since matrix diagonalization (and inversion) grows with power 3 in the dimen- 
sion, the technique of iterated resolvents might thus even be faster. - The iterated 
and the Tamm-Dancoff resolvents are distinctly different in the following aspect: 
Gn conserves particle number, but Gq does not. Both however conserve the incom- 
ing momentum. - The notation in Eq.(|41D is very compact and will be explained 
further below. - The method applies also equally well to conventional many-body 
problems since a pair-interaction generates for example the bloc matrix structure 
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shown in Figure ^ - Finally, it should be emphasized that the higher sector wave- 
functions can be retrieved by matrix multiplications from the eigenfunction in the 
lowest sector, from (l|^l/). No additional matrix diagonalizations or inversions are 
required. To show this, consider Eq.p^) for n = 1, i.e. 



{2\^) = {2\G2H2\1){1\^). 

Both 6*2 and H2 were calculated as a matrix for getting down to the effective 
interaction in the 1-space. Required is thus one additional matrix multiplication. 

Next get (3|^) = {SlG^H^ll) (1|*) + {3\G3Hs\2) (2|^). Substituting (2|^) gives 

(3|^) = (3|G3i^3(l + G2i/2)|l)(l|^). 

The general case, 

{n\^) = {n\GnHr,{l + Gn-iH,-i) . . . (1 + G^H^m (1|^), (45) 

can be proven by induction. 

6.3 A simple numerical example 

It might be interesting to study a Hamiltonian with a tridigonal band structure, 
such as was given for example in Table |: 



f{l\l + S\l) (1|F|2) . . ^^ 

(2|F|1) (2|1 + 5|2) (2|F|3) 

(3|F|2) (3|1 + 5|3) (3|F|4) 

V ■ • (4|F|3) (4|l + 5|4)y 



(46) 



According to the rules it develops the structure a of continued fraction, since with 
Gn{uj) = l/{uj — Hn) one has exphcitly 



H^ =(4|1 + 5|4), (4|^)= (4|G4FG3FG'2F|1) (1|^), 

H^ =(3|l + 5 + FG4F|3), (3|^)= (3|G3FG'2F|1) (1|^), 

H2 =(2|l + 5 + FG3F|2), (2|^)= (2|G2F|1)(1|^), 

Hi ={l\l + S + RG2R\l), (1|^). 



(47) 



Here and above a very compact notation is used, which is explained by the example 

(3|D + VG^V\'i) = (3|D|3) + (3|F|4) (416-414) (4|F|3). 

By reasons of space, the abbrevation {i\D\i) = {i\T + S\i) is used sometimes in 
the diagonal blocks. The kinetic energies T are diagonal matrix elements, and the 
seagulls do not change particle number by definition. Even that the above is very 
compact. Here is what it means in terms of matrix operations: 

{3,i\D + FG4F\3,i') = {3,i\D\3,i') 

+ Y.{^,t\F\4,j) {4,j\G,\4,k) {4,k\F\3,t'). 
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Figure 8: The the energy func- 
tion E{uj) is plotted versus uo. 
The dotted hne represents E = u. 



The 4x4 matrix 



1 



1 
2 
3 



V- 



3 

4 
5 



5 
6/ 



has eigenvalues Ei = -1.872, -0,01518, 3.333, 
and 10.55. Rs energy function is 



E{u) 



1- 1 



3-3 



UJ 



5-5 



u 



LU 



6 



The solutions of E{u) = u agree with Ei to 
within computer accuracy. The are marked 
with • in the figure. 



Within each sector \i), the number of basis or Fock states is finite, i.e. \i) = \i'.,j) 
with j = 1,2,..., Ni. The steps to be taken are then: Take the reactangular 
block matrix (4, k\E\l,i'); matrix-muhiply it from the left with the inverse matrix 
of (4|if|4) which is (4, j|G4|4, /c); matrix-multiply the result with the transpose 
(and complex conjugated) matrix {4,j\F~^\l,i)] the result is the quadratic matrix 
{3,i\FG4F\3,i')] add this to the quadratic matrix {3,i\D\3,i'). As a net result, 
the old matrix {3,i\D\3,i') in Eq.(^8D is replaced the effective (and cj-dependent) 
matrix {l,i\D{uj)\l,i'). 

To play the game a little further, one can reduce the number of states within a 
block to 1. Eq.(^6|) stands then for a tri-diagonal matrix and Eq.(^7|) stands literally 
for a continued fraction. Such one is analyzed in Figure p, with a perfect agreement 
between the method of iterated resolvents and a direct numerical diagonalization. 



6.4 The 4 by 4 block matrix for gauge theory 

In the front form, the Fock-space for harmonic resolution K = 2 has only 4 Fock- 
space sectors, see Table ^. K = 2 provides thus a good example for all 4 by 4 
hermitean bloc matrices, in which the block matrix element (2|iJ|4) vanishes. The 
case is considered as an exercise and a non-trivial example. 
Since N = 4 one starts out with the block matrix 



(4^ 





1 = (qq) 


2 = (gg) 


3 = (qqg) 


4 = (qqqq) 


l = {qq) /{1\D\1) 


{1S\2) 


(1|V3) 


(1|F4) \ 


2 = [gg) 


{2\S1) 


{2D 2) 


(21^3) 





3 = (qqg) 


(3 1-11) 


(31^2) 


{3D 3) 


(31^4) 


4 = (qqqq) 


\{AF\1) 





{MV3) 


{ADA) J 
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The block matrix element (2|iJLc|4) is a rectangular zero- matrix, since the light- 
cone operator has no matrix elements between two gluons gg and two gg-pairs. 
The notation keeps track of the field theoretic property, that each block matrix 
has only one type of interactions, namely instaneneous seagull (S) and fork (F), or 
dynamic vertex (V) interactions. 

The Hamiltonian in the 4-sector is quadratic and has the resolvent 

G4H = ,., \.., , {MH,\A) = (4|T + S\A). (49) 

Using Eq. (HI) to reduce the block matrix dimension from 4 to 3 gives 



1 = (gg) 2 = {gg) 3 = {qqg) 

1 = iqq) f {l\D + FG,F\1) {l\S\2) (1|(1 + ^^4)^13) ' 

2 = igg) {2\S\l) {2\D\2) {2\V\3) 

3 = (qqg) V (3|l^(l + G4F)|1) {3\V\2) {3\D + VG^VlS) 

Almost every block matrix element is replaced by an 'effective' element. 
Now continue to reduce from 3 to 2: The resolvent of (3|-ff|3) is now 

^3(0;) = ^3j^_^^^j3^ , (3|i/3|3) = (3|r + ^ + VG,V\3). 

With {1\A\1) = {1\T + S + FGiF + (1 + FG^)VG3V{1 + GiF)\l) one gets 

1 = (qq) 2 = (gg) 

1 = (qq) ( {1\A\1) {l\S + (1 + FG^)VG^V\2) \ 

2 = {gg) \ {2\S + ^^31^(1 + ^4^)11) {2\T + S + VG:,V\2) ) 

Finally, evaluate the last resolvent 

^2M=A^Arm' H, = {2\T + S + VG,V\2), 
to end up with the 1 by 1 block matrix 

(l|/7i|l) = {l\D^FGiF + {l + FGi)VG:iV{l + GiF)\l) (50) 

+ (1|(S + VG-iV + FG^/G^V) G2 {S + VG^V + ^6-3^04^)11). 

Here the procedure stops. For gauge theory one can order the result with the power 
of the coupling constant 



(lli^ill) = {\\T ^VG^ ^VG^VG2VG-^\1) 

+ l,\\SG2VG^VG^F ^FG^VG^VG^S ^- 

VGiVG2VG^VG^F + FG^VG^VG2VG^V + 
FG^G^VG^F + FG^VG^VG^VG^VG^F 1 1) . (51) 
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Figure 9: A particular time- ordering of the two-gluon- 
annihilation interaction is drawn in the upper left of the 
figure. If all intrinsic lines are understood as a sum 
of a 'dynamical' and an ' instantaneous' line, one gets 
the seven diagrams in the figure, drawn in two different 
conventions. One observes that the very same strings 



appear in VGSVG2VG3V as defined by Eq.(|5^. 



"1HM\^.. 



Figure 10: Each propagator Gn respresents an infinite 
resummation of Tamm-Dancoff perturbative diagrams. 
The graph (llVG^Vll) respresents a resummation of all 
Tamm-Dancoff graphs with at least one qq pair, and 
where a first gluons is emitted by the quark and a (per- 
haps other) gluon is arbsorbed by the anti-quark. 



with the abbreviations 



VG3V = VG3V + S, 
VG^VGiVG^V = VG3VG2VGsV + SG2VGsV + VG3VG2S 

+ FG^VGsV + VG3VG4F + SG2S + FG^F. (52) 



Finally, for completeness, the higher sector wave function in Eq.(|i3|) are written 
out explicitely 



(2|^) 
(3|^) 
(4|*) 



G2 

Gr4 



VGsV + VGsVG^F |1) (1|^) 



V + VGoVG-,V + VGoVG^VG.F 



|1)(1|*), 



(53) 
(54) 



VG3V + VGsVG2VG3V + VGsVG2VG3VG4F\ |1) (1|^), (55) 
with the additional abbreviations 



VG^V 



VG2VG3V 
VG3VG2VGW 



VG3V + F 

VG2VG3V + VGiF + VG2S 

VG3VG2VG3V + VG3VG4F + VG3VG2S. (56) 



These abbreviations are not very relevant for the general case. But for gauge theory 
they have a deeper physical meaning as to be discussed further below. 
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6.5 Discussion and gauge invariant interactions 

All effective interactions have a finite number of finite strings like V GV . . . GV. 
This is a colossal simplification as compared infinite number of terms of Tamm- 
Dancoff perturbation theory. Iterated resolvents resume them in a systematic way. 
Both kinds of series can be identified term by term if one expands the iterated 
resolvents also about the kinetic energies. Denoting the free resolvent by G„ it is 
related to the Tamm-Dancoff resolvent G by 

If one then expands G„ using Hn = T^ + Un one gets 

thus infinite series of infinite series, which after reordering become identical with 
the Tamm-Dancoff series. 

A look back to Eq. (|5ll) is helpful. The various term have been arranged there 
by order of the coupling constant. In the first line one finds the combination 



VG^VG2VGiV , defined in Eg. ([5^). In the caption to Figure || it is explained how 
one can simplify the formal procedures, if one associates with each intrisic line a 
dynamic plus an instantaneous line. Doing that replaces the above sum of six terms 
into a single one namely V G^V G2V G'^V , and similarily one replaces the term in 
the first line by VG^^V + S^ VG^V. 

There is a physical reason behind that: When calculating scattering diagrams 
to a fixed order in perturbation theory, one obtains the gauge invariant Feynman 
amplitudes only, if all instantaneous graphs of the same order are included. An 
example for that was given above, when calculating the scattering graph up to 
second order. Obviously, the method of iterated resolvents accounts for this aspect 
automatically, in a formal way. 

Since the remaining terms in Eq.f^TD are of order 6 or 8 in the coupling constant, 
one concludes that the restriction to the phase space for i^ = 2 satisfies gauge 
invariance only up to terms of order 6 in the coupling constant. One can verify, 
that the order of violating gauge invariance is pushed up if one includes more Fock- 
space sectors. This rule was checked explicitly in many rather laborious calculations 

n. 



m 



Theorem 1 An easy trick to achieve gauge invariance is the following: Set to zero 
formally all instantaneous interactions, perform the steps required by the resolvents. 
At the end, when having calculated the effective interaction in terms of the vertex 
interaction V , one replaces each internal line by the sum of a dynamical and a 
instantaneous line, according to the general rules of light- cone perturbation theory. 
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If one follows tlies rules for the 4 by 4 matrix given in Eq. (|48|) one obtains for 
the sector Hamiltonians and wavefunctions simply 

H4 = T4, {A\^) = {AlG^VGaV + G4VG3VG2VG3V\1) {1\^), 

Hs = n + VG^V, {3\^) = {3\G3V + G3VG2VGsV\l){l\^), 

H2 = T2 + VGsV, (2|*) = (2|\/G'3V^|1)(1|^), 

H^ = T1 + VG3V + VG3VG2VG3V. (57) 

These expressions are more transparent than the full equations; they are also gauge 
invariant. 

6.6 The iterated resolvents for arbitrary K 

When one substitutes formally all instantaneous interactions by block matrices with 
only zeros, keeping only the vertex interactions V, one gets a block matrix structure 
as displayed in Table ^ Because of the many zero matrices, the construction of 
the the sector Hamiltonians is now rather straightforward and simple. Here they 
are for the first 12 sectors: 

H12 = T^2 + VGnV + VGnVGiGVGi2V + VGi3V, (58) 

Hu = Tn + VGieV + VGiQVGi5VG,eV + VG^2V, (59) 

Hio = Tio + VGi5V + VGi5VGuVGi5V + VGiiV, (60) 

H, = T, + VG,oV + VGuV, (61) 

Hs = T8 + VGi2V + VGi2VGnVGuV, (62) 

Hr = Tj + VGnV + VGuVGioVGuV + VGsV, (63) 

Hq = Tq + VGioV + VGioVGgVGioV + VGjV, (64) 

H5 = T, + VGeV + VGgV, (65) 

H4 = Ti + VGrV + VG7VG6VG7V, (66) 

H3 = Ts + VGeV + VGeVG5VGeV + VGiV, (67) 

H2 = T2 + VGSV + VG5V, (68) 

H^ = T1 + VG3V + VG3VG2VG3V, (69) 

Actually, for K = 4, the exact expressions corresponding to Table | are obtained 
by setting 6*13 = l/{uj — T13) and Gn = for ra > 14. But it is easy to write down 
the expressions in Eqs.(pB[)-(|5DD valid for a higher value of K. One also notes that 
they give all sector Hamiltonians for ii' = 2 in Eqs.(|57|), by setting formally Gn = 
for n > 5. 

Theorem 2 For sufficiently large harmonic resolution K the formal expressions 
for effective Hamiltonian in sufficiently low sectors become independent of K. 

One therefore can go to the limit K ^ 00 and thus to the continuum limit. 
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Table 8: The Fock space and the Hamiltonian matrix H' = T + V for a meson at 
Gxed value ofK = 4. — See discussion in the text. The diagonal blocks are denoted 
by T. Most of the block matrices are zero matrices, marked by a dot (■). The block 
matrices marked by V are potentiall non-zero due to the vertex interaction. 



6.7 Propagation in medium 

Here seems to be a problem: For calculating G3 one needs Gq, G5 and G4, for 
calculating Gq one needs Gio, Gg and G7, and so on. In the next section will be 
shown how the hierarchy can be broken in a rather effective way. That final step 
will be comparatively simple if one has analyzed the propagators for the sectors 
with one gg-pair and arbitrarily many gluons, as follows next. 

Consider first the the effective interaction in the space of one gg-pair and one 
gluon as given by Eg . (|67D . The corresponding diagrams can be grouped into two 



topologically distinct classes, displayed in Figs. [T^ and |T^. The the diagrams in 
Fig. |l^ are obtained by adding a non-interacting gluon line to the former diagrams 
in Fig. |rT|. The gluon does not change quantum numbers under impact of the 
interaction and acts like a spectator. Therefore, the graphs in Fig. will be 
refered to as the 'spectator interaction' U3. In the graphs of Fig. |l3]the gluons are 



scattered by the interaction, and correspondingly these graphs will be refered to as 
the 'participant interaction' U3. Thus, U3 = U3 + U3. The separation into a graph 
with only one interacting quark-pair Fock-space sectors 



Un = Un + Un 



(70) 



except in those with only gluons. More explicitly, the spectator and participant 
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Figure 11: The three graphs of the ef- 
fective interaction in the qq-space. The 
vertical hnes denote the propagators 
Gn- On the right, the vertical lines de- 
note the propagators Gn- The coupling 
function at the vertices is symbolized by 
graphs as they would appear in a per- 
turbative analysis. 





interactions in the lowest sectors with one quark-pair become 

U3 = VGqV + VGeVG^VGeV, f/3 = VG^V + VGqV, 



Ue 



VG.nV + VGiaVGgVGioV, 



nov, 



f/fi 



VG7V + VGioV. 



(71) 



The same operators can appear in both U and U, but they refer to different graphs. 
Since the Hamiltonian is additive in spectator and participant interactions, t/„ can 
be associated with its own resolvent 



Gn 



( while Gn = 



U-Tn-Un V UJ-Hn CU - T„ - f/„ - f/„ 

The relation of Gn to the full resolvent is Gn = Gn + Gn Un Gn, or 

Gn = Gn + GnUn Gn + GnUn GnUn Gn + • • • , 



(72) 



(73) 



writing it as an infinite series. Note that the propagator Gn contains the interaction 
Un in contrast to the free Tamm-Dancoff propagator G in Eg. (PD . One deals here 
therefore with 'perturbation theory in medium'. But since it is driven to all orders 
one better speaks of 'propagation in medium'. 
The series can be identically written as 



Gr 



"-^n ^^n '^n, 



with Gn 



1 — GnUn 



(74) 
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Figure 12: Tie tiree graphs of the Figure 13: Some six graphs of the par- 
spectator interaction in the qqg-space. ticipant interaction in the qqg-space. 
Note the role of the gluon as a spectator. 

One can verify this order by order in perturbation theory or by the identity 



{U - Hn){l - GnUn) = {UJ-Hn)[l ^Ur, = CU - f^n- (75) 

The inverse gives C\Gn = Gn, and with the identity GnGn = GnC^, one gets 
Eg. (|7iD . This remarkable property is pecuhar to the method of iterated resol- 
vents. Whenever a quark-pair-glue resolvent is sandwiched in between two vertex 
interactions, it allows to define a modified vertex interaction V, since 



VGnV^ = VGn Gn C^ V^ = VGn V\ with V = VGn- 

One can thus rewrite the sector Hamiltonians like for example 



f/e = VGioV + VGioVG<,VGioV, 
Us = VGeV + VGeVG^VGeV, 
U^ = VG3V + VG3VG2VG3V. 



(76) 



(77) 
(78) 
(79) 



Now they have all essentially the same structure, contrary to Eqs.(p5D where they 
were different. Note that G2, G5 and Gg are pure gluon propagators. 

Generally spoken, the rectangular block matrix V is multiplied with the square 
matrix C„. Below it will be shown that C„ is approximately diagonal and indepen- 
dent of the spin. Each vertex matrix element is multiplied with a simple function 
which depends on the momentum transfer Q across the vertex. Equivalently one 
replaces the coupling constant g hj g = gCn{Q), such that C„ can be interpreted 
as a coupling function. 
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6.8 The exact effective interaction 

The most important result of this section is that gauge theory particularly QCD 
has only two structurally different contributions to the effective interaction in the 
gg-space, see Eq.(^). It scatters a quark with helicity \g and four-momentum 
p = {xP^,xP± + k_j_,p~) into a state with A' and four- momentum p' = [x' P~^ , x' P± + 

k^,p'~), and correspondingly the antiquark. In the continuum limit, the resolvents 
are replaced by propagators and the matrix eigenvalue problem Hes\ip) = M'^ltp) 
becomes an integral equation with a rather transparent structure: 



M^{x,kr,Xg,Xg\tpi) 



nig + kl m^ + kl 
X 1 — X 



(a;,4;A„A,-|^i) (80) 



+ Y. fdx'd^k[ R{x\ fci) {x, kr, A„ Xg\UoGE\x\ k[; a;, a;.) {x', k[; a;, A'^I^,) 

+ J2 dx'd'^k^ R{x', k'^) (x, k^\ \, WU-^gaW , k[; A^, A'^) {x, k'j_; A^, X'glipi). 



K'K 



The domain of integration is set by the cut-off function R given in Eq.(]35[). The 
effective one-gluon exchange 

f/oGE = VGsV (81) 



conserves the flavor along the quark line. As illustrated in Fig. |1T] the vertex 
interaction V creates a gluon and scatters the system virtually into the qqg-spa.ce. 
As indicated in the figure by the vertical line with subscript '3', the three particles 
propagate there under impact of the full Hamiltonian before the gluon is absorbed. 
The gluon can be absorbed either by the antiquark or by the quark. If it is absorbed 
by the quark, it contributes to the effective quark mass m. The second term in 
Eq.(|79D is the effective two-gluon-annihilation interaction, 



f/TGA = VGsVG2VGsV, (82) 



as represented by the graph Uifi in Fig. 11. The virtual annihilation of the qq- 



pair into two gluons can generate an interaction between different quark flavors. - 
The eigenvalue M^ is one of the invaraint-mass squared eigenvalues of the full 
light-cone Hamiltonian, and the wavefunction {x, fcj_; A^, Ag|^j) gives the probability 
amplitudes for finding in the gg-state a flavored quark with momentum fraction x, 
intrinsic transverse momentum k_j_ and helicity Ag, and correspondingly an anti- 
quark with 1 — a;, —k± and Xg. 

7 The breaking of the propagator hierarchy 

The content of the preceeding sections is exact but rather formal. In the sequel, 
rigor will be given up by the aim to obtain a solvable equation. In order to spot 
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easier the essential approximations, they will be dressed as 'theorems', which can 
- or cannot - be proven later. 

All sector Hamiltonians can be diagonalized on their own merit, for example 

Ml{q;q\<il,.,) = Y.i9;q\Hi\q';q) {q;q\^i;), 
Ml{q;q;g\^,,,) = E (?! ^S ^l^slg; g; ^') (?';?; ^1^3;.). (83) 

The spectra might be continuous, but the eigenvalues are denumerated for sim- 
plicity by i. By construction, one knows the relation between these two sets of 
eigenvalues, by the following reason: Because of the separation into spectators 
and participants the gluon in the qq g-equation is a non-interacting particle, it 
moves freely relative to a gg-bound state. If its four-momentum is parametrized as 
q'^ = {yP^, yP± + qL,qg) one has 

Ml, . £!±| , C. (84) 

Every qq bound state M is a band-head in the qqg spectrum. With uj = MfQ{uj) = 
M^ one gets 

. _ Ml, . M' - ('4^ . «) ^ -«^^f±^. (85) 

''' \ i^-y) y) y{i-y) ^ ^ 

Knowing the eigenvalues and the eigenfunctions of the qq bound state one alos 
know the spectrum and the eigenfunctions of H^. Knowing the eigenfunctions, the 
exact resolvent can be calculated. 

Theorem 3 The exact propagators can be approximated by closure. 

The substitution of the exact propagators by closure is a widely used approximation 
in many-body theory and in chemistry. In principle, the exact propagator is non- 
diagonal. Performing closure, it becomes diagonal, which is reflected by the Dirac- 
delta function (g; q; g\q'; q'; g') in the single particle momenta and helicities 

(g; g; ^l^slg'; g ; g) = (g; g; g\q'; q; g')G3{q; g; g). (86) 

The Dirac-delta function is multiplied with the function 



y(i - y) 

y^M^ + ql 



G3{q;q;9) = - Jl,, L - ^ 



The in-medium propagator G^ ceases to be a functional of the higher resolvents. 
The hierarchy of the iterated resolvents is broken. 
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Figure 14: The qq vacuum polarization Figure 15: The gg vacuum polarization 
graph. graph. 

Theorem 4 In the solution, the interacting particles propagate free particles. The 
in-medium propagators can be replaced by the free propagators. 

The free propagator in the qqg-space can be written 

1 y y /ooN 



a 



3, free 



P+{p- 



P 



Q' 



Q' 



t/2(2m)2 + ql 



•2' 



where Q"^ = {p — p'Y is the four- momentum transfer along the quark hne, in the 
notation of Fig. |14[ For sufficiently small y holds {p—p'Y = ~[i/^(2^)^ + 5l]. If one 
substitutes M ~ 2m, which holds to rather good approximation, the momentum 
transfer in Eg . (^Sf ) is the same as in Eq.(|87|). 

Similar considerations also hold all other propagators with at least one gg-pair. 

Having replaced the in-medium-propagators Gn by the free propagators, one 
can calculate the graphs like in the usual (light-cone-) time-ordered perturbation 
theory. The diagram Ui 2 in Fig. O yields the effective quark mass 



— 2 I I 2 



, a n? — 1 , A 



71 2nr 



In— r. 



m„ 



^9) 



Correspondingly one can calculate the coupling function 



C, 



^ ^ ~ = 1 + ^Gst/s + ^GsUsGsUs + ... 

1 — G3U3 ^ o 



(90) 



as defined by Eq.(|74D. With f/3 = VG4V + VGqV from Eq.(|7lD, and thus U^ 
VG4V + VGqV from Eq.(^), the first two terms in the expansion are 



Vc^V+ -{VG3VG4V + VG^VGeV) ~ V + -(^^31^^41^ + VGaVGeV). (91) 

In the last step R4 = Rq = 1 was set. In general, they contribute terms of higher 
orders in g, which must be suppressed for consistency. Figs. |l^ and ^ give only 2 
representative graphs. According to the rules of time ordered perturbation theory, 
there are about 22 different time-ordered and instantaneous diagrams. All of them 
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were re-calculated by Raufeisen [Q, for all values of Q. As shown to some detail 
in 



33| one gets for sufficiently large A 



^ TT/^ ^ ^Tr llaric , ( A^ \ a / A^ 

G,V{G, + G,)V = ^T^ln — ^-7^ -^Eln 



127r \Aml + Q^J Gtt y Kim} + Q^ 
= a6oln(AV«:2)-6(Q), &o = ^^^^V^' (92) 

iZTT 

with an arbitrary k, and a ^(Q) independent of A given in Eg . (|93|) . One conculdes 
that the coupling function C^ is identical with the familiar vertex correction. In the 
limit where Q is larger than all mass scales, this result agrees with the calculations 
of both, Thorn |3^ and Perry |3^], having been done previously. 



What happens if one substitutes the free propagators in Eq.(pO|) by the non- 
perturbative propagators G4 = (cj — Hqqqq)~^ and Gq = {uj — Hqqgg)~^, at least in 
an approximate fashion? — There are additional graphs. In the fermion loop of 



vacuum polarization appear two graphs in addition to Fig. |lj. In one of them, a 
gluon is emitted and absorbed on the same quark line which changes the bare quark 
mass TUf into the physical quark mass rrif. In the other graph, the gluon is emitted 
from the quark and absorbed by the anti-quark which represents an interaction. 
In consequence one has a bound state with a physical mass scale fif. We replace 
therefore 2m f =^ 2m/ =^ fif. Similar considerations hold for the gluon loop in 
Fig. |15| and lead to the substitution 2mg =^ 2m g =^ fig. Both fig and fif are 
interpreted as physical mass scales. The physical gluon mass nig vanishes of course 
due to gauge invariance. This is not in conflict with f.e. Cornwall's suggestion of 
a flnite effective gluon mass |^ since one can deflne (mg)eff = /ig/2. 



As a consequence of the approximation in Eq. (|86|) , the propagators G3 and G4 
are diagonal in Fock space, which in turn leads to diagonal products G3VG4V. 
They can be resumed therefore to all orders according to Eq.(^). With a = 
g^Cl/A'K, one gets 

a(Q;A) = — , , , ,0 1 o\ 777^ with 

The effective flne structure constant depends on the momentum transfer Q across 
the vertex and on the cut-off A. Now, all the pieces are together which are needed 
for a further discussion of the effective interaction as deflned in Eq.(^^). Finally, 
the the instantaneous interaction in the effective one-gluon-exchange interaction 
UoGE = VG^V is restored according to VG-iV — ^ Cl{VG:iV + S), with the 
expression given in Eq. (p6|) . 
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8 Status and discussion 

What have we reached? - The general expression for the effective Hamiltonian in 
the gg-sector was given in Eg. (|80D . If one restricts to consider mesons in which the 
quark and the anti-quark have different flavors, fq ^ fq, the two-gluon annihilation 
diagram can not get active and one remains with 



M^{x,k_i_;Xq,Xq\i'n) 



m2(A) + kl m?(A) + kl 



+ 

X i — X 



{x,k_i_;Xq,Xq\lpn) 



+ J2 fdx'd^k'^ {x, kr, Xq, Xq\UoGE\x', fci; A;, A',) {x', k'^- A;, A'.lV'n) • (94) 

The essential achievement of the two preceeding sections is that the kernel of this 
integral equation can be written down in the very explicit form 






(Ag, Ag|f/oGE|Aq, X'q) o^2n2 ^"^9' Xq\S\\, \) ^ ^^ \ ]^^ =■ (95) 



All many-body effects reside in the effective mass mq(A) and in the effective cou- 
pling constant a(Q;A). They are given in Eqs. (|89|) and (|93|) , respectively. The 
mean of the four-momentum transfers of quark and anti-quark, Qi = —{kg — fc')^ 

and Qg = {kg — k'^f, respectively, is denoted by Q^. The spinor factor is 

(a„a|5|a;,a') = [u{k,\)ru{k\>:)\ [u{k,\)^,u{k\\%. 

The cut-off function R sets the domain of integration and is defined in Eq. (|35D . 
The approximations made have been carefully enumerated in Section 7. 

One has thus reached the goal of reducing the field theoretical many-body 
problem of M^ = Hi^q to a one-body problem with an effective interaction. For 
any value of the Lagrangian coupling constant g, the flavor masses m/, and the cut- 
off A, one can generate the eigenvalues M^ on a computer, very much in analogy 



to QED [^, ^. With the known eigenfunction in the gg-sector one an generate 
all higher Fock-space amplitudes, according to the explicit formulas in Section 6. 
It looks as if one has solved the problem. 

But one has solved the problem only in a superficial way, since the solutions 
depend on the unphysical parameter A. One of the central issues in gauge-field 
theory is to remove this dependence by a renormalization group analysis, see below. 

The eigenvalues of the integral equation are the invariant masses squared. Very 
often it is more convenient to think in terms of an analogue of a non-relativistic 
Hamiltonian Hs-, 

Hlc = {mq + mqiY + 2{mq + rriq') Hs, (96) 

whose eigenvalues E = Hs are the more convential binding energies, while the 
eigenfunctions are the same. As shown in |0, the so defined 'Schrodinger operator' 
has much in common with a non-relativistic Hamiltonian for the Coulomb problem. 
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8.1 Renormalization group analysis 

All eigenvalues M^ depend thus on the cut-off A, the bare masses m and the bare a. 
In line with modern interpretation of quantum field theory [|^ , they are unphysical 
parameters and can be replaced by (mA,aA), by functions of A, such that the 
eigenvalues are cut-off independent, i.e. 

^M2(A,«A,mA)=0. (97) 

This fundamental equation of the renormalization group must hold for all eigen- 
values. 

The effective Hamiltonian H = Hi^ceS depends on A only through rfif, a and 
the regulator function R. The renormalization group equations can be written 
therefore in the operator form 

Smj ^m + fe^l^n) + 5i?^|^n) = 0. (98) 

omj da on 

The simultaneous variation of all three terms is a very difficult problem. But one 
can replace this equation by three independent ones, 

^ruf = 0, (99) 

-^a = 0, (100) 

^jj^m = 0. (101) 

This is possible since one still solves Eq.(^Sj). 

Eqs.(|99D and (|100|) are two equations for the two unknown functions ttt-a and 
a\. The first equation says that the effective flavor masses rrif = rn/(A, oa, "^a) 
are renormalization group invariants. The numerical value of mj must be fixed by 
experiment. - Eq.( |100| ) is then considered as an equation for a\ at fixed values of 
rrif. Using the expression for a{Q; A) given in Eq.(p3D yields then |^ 



"^ = 6^m1v^' '^"^ ''^^^ = b(Qy ^'°'^ 

with b{Q) given in Eg . (|93D . All A-dependence cancels exactly in favor of the renor- 
malization group invariant k, which is sometimes called the QCD-scale Aqcd- The 
scale K must be determind from experiment. - Once a^ is known, Eq.(^) is an 
equation for ttia- This function can be determined from Eq. (PU]) , but we renounce 
to do that here, since m\ is needed nowhere in the present context. - Having 
fixed the functions oa and ttia one has exhausted all freedom provided by conven- 
tional renormalization anaysis of field theory [|^ . Since the cutoff-dependence was 
removed by renormalization, the cut-off can be driven to infinity, thus 

R{x,kr,A) = 1. 
Formally, this solves Eq.( p.01| ). 
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8.2 The remaining divergence 

After renormalization, the kernel of the one-body equation is independent of A. 
Let us discuss its relevant part 

1 S 

At very small momentum transfers Q^ -^ 0, where x ~ x' and kj_ ^ k'j_, the effective 
coupling constant locks into the finite value b{0), see Eg . (|93D . The spinor function 
also goes to a constant, S ~ Afngfrig, see |^T| or 0. The remaining 'Coulomb 
singularity' Q"^ is square-integrable and can be dealt with by convenient numerical 
means ||3^, |3l|]. For very large momentum transfers Q^ -^ oo, the behaviour is 
quite different. Both S and Q^ tend to diverge, S oc (/^l)^ and Q^ oc (fc|)^, 
but such that the ratio S/Q"^ tends to a finite constant. Disregarding the very 
slow logarythmic increase of b{Q), the kernel of the integral equation is therefore 
essentially a dimensionless constant 



/C ~ constans, for Q^ ^ 0. 



One has sufficient evidence [31] that this behaviour creates all kinds of problems, 
among them a diverging eigenvalue. 

One can separate the kernel of the integral equation identically into two pieces 
/C = /Ci + /C2, i-G. 

IS fi^ + Q^ ^ l S S 



2V 



b{Q)Q^ IJ^ + Q^ 6(g)g2(i + g2/^2) 5(g)(^2 + g 

with an arbitrary mass parameter fi. The first part /Ci is well behaved and vanishes 
at least like Q~^ in the limit of large Q. The second part /C2 is a strange object: it 
is a constant almost everywhere, in particular 



S_ 



K.2 ~ constans, since — ■ ~ constans, for Q^ ^ 0. 



One should emphasize the following aspect. The typical field theoretical diver- 
gences like the divergence of the effective coupling constant lead to a divergence 
of the kernel, but adding a finite constant to the kernel leads also to a divergence. 
Thus far one does not understand the reason for the latter. 

Perhaps one can get further insight by studying the following eigenvalue in 
(usual) three-momentum space (/c) 



-f k -> cv r 
E^(k) = —i:{k)--d'k' 
2m, TT"^ J 



1 1 

{k - k'Y fi^ 



ip{k'). (103) 



In the spirit of Eq. (PU]) , the right-hand side is a conventional Schrodinger equation 
for the Coulomb problem with the reduced mass 1/m = l/mg + l/mg, but where an 
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Table 9: The empirical masses oi the Havor-off-diagonal physical mesons in MeV. 
Vector mesons are given in the upper, pseudo-scalar mesons in the lower triangle. 
Their physical nomenclature is given on the right. 
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Table 10: The calculated masses oi the Bavor-off-diagonal mesons in MeV, as 
obtained from a fit to Eq.l \103{ ) on the left and to Eq.( \10^ ) on the right. Vector 
mesons are given in the upper, pseudo-scalar mesons in the lower triangle. The 
* upperscript marks mesons where the shift and mass parameters have be fitted to. 



arbitrary constant has been added in the kernel. In Heidelberg, we are presently 
working on the question whether its eigenvalue E diverges. Preliminary studies 
with the Fourier transform of Eq.( |103| ), where the constant works its way into a 
weird Dirac-delta function 5^^^ (x) as a potential, 



V2 



^ + ^S''\ 



X] 



\x\ 



/x^ 



0(f), 



indicate indeed that the eigenvalue of the Is-state in particular diverges. 

If these preliminary results materialize, one must regularize and subsequently 
renormalize Eq.( |103| ). The energy of the Is-state will then be essentially a renor- 
malization constant E^ = s. Given that to be true, all of us have been riding the 
wrong horse. It is not the regular part of the interaction kernel /Ci, but it is the 
'constant' /C2 which determines the mass in the first place. Since /C2 involves a 
renormalization and since the singlet can have an energy shift (s_) different from 
the triplet (s+), one has at most two additional renormalization constants, subject 
to be determined by experiment. 
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Table 11: The quark mass parameters m,q and the energy shifts s±, all in MeV, 
as obtained from a fit to the physical meson masses. 

8.3 A mass formula 



Is there some physical evidence for the above considerations? - Inspired by Eq. (|96D , 
one can seek the masses of the physical mesons in the form 13^ 



M' 



{mq + mqiY + 2(mg + mg/) s± 



(104) 



and fit the energy shifts s± and the physical quark masses fUg to the physical meson 
masses, compiled in Table ^ from the data of the particle data group [^ which do 
not include yet the topped mesons. 

The following procedure was applied. First, the up and the down mass were 
chosen equal. Then the empirical masses of vr'*' and p+ were used to determine the 
mass shifts for the singlet and the triplet. The remaining quark masses are obtained 
from the pseudo-scalar mesons with an up quark, which exhausts all freedom in 
determining physical parameters. The remaining 13 off-diagonal pseudo-scalar and 



vector meson masses are calculated straightforwardly and compiled in Table 10 



In comparing the so obtained numbers with the experiment, one notes with great 
surprise that the discrepancy exceeds only rarely an estimated error of 10%. The 
resulting quark masses are given in Table |ll|. 

Actually, there is no reason why one should stick to Eq. (|104|) . Equivalently, one 
can replace s± by s± a/ijriq + nig'). Choosing a = m^ + m^ = 700 MeV, one gets 



M 



2 l^^-r 



(nig -h mgt) -h 2(ms + mi) s±. 



(105) 



A form like this was suggested to me by Dae Sung Hwang during the lectures. 
Redoing the fit for the heavy quarks gives the results compiled in Table |10|. The 



discrepancy with the experimental numbers is now on the level of 1%. This excellent 
agreement should not be over-emphasized. In any case it does not falsify the above 
considerations. 



9 A short summary 

In these lecture notes the attempt was made to phrase and discuss a relativistic 
quantum theory such as QCD from the point of view of Dirac's front form of 
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hamiltonian dynamics. Sometimes refered to as light-cone quantization, it has the 
goal to describe and understand the bound-state structure of hadrons in terms of 
their constituents from a covariant theory. As the lectures show, this goal has not 
been reached yet. 

As reviewed in [^, there might be alternative roads to reach the same goal, 
like for example the Bethe-Salpeter equations, lattice gauge theory, Hamiltonian 
flow equations, just to name a few. Each of these approaches have their own 
inherent adventages or disadventages. By reasons of space their discussion is left 
out in these notes. By the same reason the light-front renormalizations approach of 



Wilson and collaborators |^, ^ is omitted. The many articles have thus far not 
provided conclusive evidence for success, apart from the fact that virtually every 
Langrangian symmetry has been violated in the approach. It is much too early to 
draw definitive conclusions. 

The front form is however useful to study the structure of arbitrary field the- 
ories. The theories considered are often not physical, but are selected to help in 
the understanding of a particular non-perturbative phenomenon. The relatively 
simple vacuum properties of light-front field theories underly many of these 'an- 
alytical' approaches. The relative simplicity of the light-cone vacuum provides a 
firm starting point to attack many non perturbative issues. As we have seen, the 
problems in two dimensions, in one space and one time dimension, not only are 
they tractable from the outset but in many cases like the Schwinger model can they 
be solved numerically by 'discretized light-cone quantization (DLCQ)' to almost 
arbitrary accuracy. This solution gives a unique insight and understanding. The 
Schwinger particle indeed has the simple parton structure that one hopes to see in 
QCD. 

Unfortunately the same cannot be said for the physical world in 3 space and 1 
time dimensions. The essential problem is that the number of degrees of freedom 
needed to specify each Fock state even in a discrete basis grows much too quickly. 
As discussed in this review, the basic procedure is to diagonalize the full light-cone 
Hamiltonian in the free light-cone Hamiltonian basis. The eigenvalues are the in- 
variant mass squared of the discrete and continuum eigenstates of the spectrum. 
The projection of the eigenstate on the free Fock basis are the light-cone wavefunc- 
tions and provide a rigorous relativistic many-body representation in terms of its 
particle degrees of freedom. Given the light-cone wavefunction one can compute 
the structure functions and distribution amplitudes. More generally, the light-cone 
wavefunctions provide the interpolation between hadron scattering amplitudes and 
the underlying parton subprocesses. The method of iterated resolvents can be a 
useful intermediary step to generate these wavefunctions. The unique property 
of light-cone quantization that makes the calculations of light cone wavefunctions 
particularly useful is that they are independent of the reference frame and that the 
same wavefunction can be use in many different problems. 

Finally, let us highlight the intrinsic advantages of light-cone field theory: 
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• The light-cone wavefunctions are independent of the momentum of the bound 
state - only relative momentum coordinates appear. 

• The vacuum state is simple and in many cases trivial. 

• Fermions and fermion derivatives are treated exactly; there is no fermion 
doubling problem. 

• The minimum number of physical degrees of freedom are used because of 
the light-cone gauge. No Gupta-Bleuler or Faddeev-Popov ghosts occur and 
unitarity is explicit. 

• The output is the full color-singlet spectrum of the theory, both bound states 
and continuum, together with their respective wavefunctions. 
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